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Abstract The quantum information approach to many-body physics has been very 
successful in giving new insight and novel numerical methods. In these lecture notes 
we take a vertical view of the subject, starting from general concepts and at each 
step delving into applications or consequences of a particular topic. We first review 
some general quantum information concepts like entanglement and entanglement 
measures, which leads us to entanglement area laws. We then continue with one 
of the most famous examples of area-law abiding states: matrix product states, and 
tensor product states in general. Of these, we choose one example (classical super- 
position states) to introduce recent developments on a novel quantum many-body 
approach: quantum kinetic Ising models. We conclude with a brief outlook of the 
field. 



1 Introduction 

There has been an explosion of interest in the interface between quantum informa- 
tion (QI) and many-body systems, in particular in the fields of condensed matter and 
ultracold atomic gases. Remarkable examples are Ref. 1 1 1, which proposed using ul- 
tracold atomic gases in optical lattices for QI (and stimulated interest in distributed 
quantum information processing), and Refs. ||2]-Q, who discussed the first connec- 
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tions between entanglement and quantum phase transitions (QPT). Overall, the con- 
fluence of ideas has opened fundamentally deep questions about QPT's, as well as 
practical questions about how to use QI ideas in numerical simulations of many- 
body quantum systems. Here, we will (partially) review these two major themes. 
We will first introduce some basic notions and tools of quantum information theory, 
focusing on entanglement and entanglement measures. We shall then discuss area 
laws, i.e. laws that characterize correlations and entanglement in physically relevant 
many-body states, and allow to make general statements about computational com- 
plexity of the corresponding Hamiltonians. Afterwards, we will explore the concept 
of matrix product states (MPS) and their generalizations (projected entangled pairs 
states, PEPS, and tensor networks states). These states provide not only a very use- 
ful ansatz for numerical applications, but also a powerful tool to understand the 
role of entanglement in the quantum many-body theory. We will review one par- 
ticular example of a state with a straightforward MPS representation: the classical 
superposition state. The introduction of its parent Hamiltonian will lead us to the 
final subject of these lectures; quantum kinetic Ising models — an analytically solv- 
able generalization of the popular classical many-body model described by a master 
equation. 



2 Aspects of Quantum Information 

Quantum theory contains elements that are radically different from our everyday 
("classical") description of Nature: a most important example are the quantum cor- 
relations present in quantum formalism. Classically, complete knowledge of a sys- 
tem implies that the sum of the information of its subsystems makes up the total 
information for the whole system. In the quantum world, this is no longer true: there 
exist states of composite systems about which we have complete information but we 
know nothing about its subsystems. We may even reach paradoxical conclusions if 
we apply a classical description to such "entangled" states — whose concept can be 
traced back to 1932 in manuscripts of E. Schrodinger. 

What we have just realized during the last two decades is that these fundamen- 
tally nonclassical states (from hereon "entangled states") can provide us with more 
than just paradoxes: They may be used to perform tasks that cannot be achieved with 
classical states. As landmarks of this transformation in our view of such nonclassi- 
cal states, we mention the spectacular discoveries of (entanglement-based) quan- 
tum cryptography |5|, quantum dense coding [6], and quantum teleportation (|7|. 
Even though our knowledge of entanglement is still far from complete, significant 
progress has been made in the recent years and very active research is currently 
underway (for a recent and very complete review see |8|). 

In the next section, we will focus on bipartite composite systems. We will define 
formally what entangled states are, present some important criteria to discriminate 
entangled states from separable ones, and show how they can be classified according 
to their capability to perform some precisely defined tasks. However, before going 
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into details, let us introduce the notation. In what follows we will be mostly con- 
cerned with bipartite scenarios, in which traditionally the main roles are played by 
two parties called AUce and Bob. Let Ha denote the Hilbert space of AUce's phys- 
ical system, and Hb that of Bob's. Our considerations will be restricted to finite- 
dimensional Hilbert spaces, so we can set Ha = C''-* and He = C''^. Thus, the 
joint physical system of AUce and Bob is described by the tensor product Hilbert 
space Hab = Ua^Hb = C^^ O C*^. Finally, B{H) will denote the set of 
boimded linear operators from the Hilbert space HtoH. 



2.1 Bipartite pure states: Schmidt decomposition 

We start our study with pure states, for which the concepts are simpler. Pure states 
are either separable or entangled states according to the following definition: 

Definition 1. Consider a pure state \tpAB) from Ha ® 'Hb. It is called separable if 
there exist pure states I ?/^^) € 'Ha and |V'b) € "Hs such that I'^ab) = IV'a) <8i IV'-Bl- 
Otherwise we say that IV'ab) is entangled. 

The most famous examples of entangled states in Hab are the maximally entangled 
states, given by 

d-l 

\i'f) = ^,y.\^A®\t)B {d = ram{dA,dB}), (D 

where the vectors and {|i)B} form bases (in particular they can be the stan- 

dard ones) in Ha and Hb, respectively. In what follows, we also use the notation 
to denote the projector onto The reason why this state is called maxi- 

mally entangled will become clear when we introduce entanglement measures. 

In pure states, the separability problem — the task of judging if a given quantum 
state is separable — is easy to handle using the concept of Schmidt decomposition: 

Theorem 1. Let \iI)ab) e Hab = C''-* ® C''^ with dA < dB- Then IV'ab) can be 
written as a Schmidt decomposition 

r 

IV-ab) =$^Ai|ei)«)|/i), (2) 

i=l 

where \ ) and \ fi ) form a part of an orthonormal basis in Ha and H b, respectively, 

Proof. A generic pure bipartite state \tpAB) can be written in the standard basis 
of Ha (8) Hb as IV'ab) = Z^i^o Si=o "^il*) ® b)' where, in general, the co- 
efficients aij form a dA x dB matrix A obeying tr(ylM) — 1. Using singular- 
value decomposition, we can write A = VDaW^, where V and W are unitary 
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(V^^^V^ = W^W = 1^) and Da is diagonal matrix consisting of the eigenvalues Xi 
of 1^1 = VAhi. Using this we rewrite IV'ab) as 

1=0 j=0 fc=l 

where r < < ds denotes the rank of A. By reshuffling terms, and defining 
|efe) = ^^^=(7^ and l/fc) = X^jfcT^ C^ifcli) we get the desired form [Eq. (B]. 

To complete the proof, we notice that due to the unitarity of V and W, vectors ]e,) 
and \fi) satisfy {ei\ej) = {fi\fj) = Sij, and constitute bases of Ha and Hb re- 
spectively. In fact, {A^, \ei)} and {Xf, \ fi)} are eigensystems of the first and second 
subsystem of IV'ab)- Moreover, since tr(yl^yl) = 1 it holds that Xf — l.U 

The numbers Xi > {i — 1, . . . ,r) are called the Schmidt coefficients, and r 
the Schmidt rank of IVmb)- One can also notice that {Af , je^)} and {A|, are 
eigensystems of the first and second subsystem of | V'ab), and that the Schmidt rank 
r denotes the rank of both subsystems. Then, comparison with definition [T] shows 
that bipartite separable states are those with Schmidt rank one. Thus, to check if a 
given pure state is separable, it suffices to check the rank r of one of its subsystems. 
If r = 1 (the corresponding subsystem is in a pure state) then \iI}ab) is separable; 
otherwise it is entangled. Notice that the maximally entangled state ([T]i is already 
written in the form Ob, with r = d and all the Schmidt coefficients equal to 1/ Vd. 



2.2 Bipartite mixed states: Separable and entangled states 

The easy-to-handle separability problem in pure states complicates considerably in 
the case of mixed states. In order to understand the distinction between separable 
and entangled mixed states — first formalized by Werner in 1989 |9l — let us con- 
sider the following state preparation procedure. Suppose that Alice and Bob are in 
distant locations and can produce and manipulate any physical system in their lab- 
oratories. Moreover, they can communicate using a classical channel (for instance 
a phone line). However, they do not have access to quantum communication chan- 
nels, i.e. they are not allowed to exchange quantum states. These two capabilities, 
i.e. local operations (LO) and classical communication (CC), are frequently referred 
to as LOCC. 

Suppose now that in each round of the preparation scheme, Alice generates with 
probability pi a random integer i (i = 1, . . . , K), which she sends to Bob. Depend- 
ing on this number, in each round Alice prepares a pure state \ei), and Bob a state 
l/i). After many rounds, the result of this preparation scheme is of the form 



K 

QAB = ^p^\ei){e^\ ® 

i=l 



(4) 
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which is the most general one that can be prepared by Alice and Bob by means of 
LOCC. In this way we arrive at the formal definition of separability in the general 
case of mixed states. 

Definition 2. We say that a mixed state qab acting on T-Lab is separable if and only 
if it can be represented as a convex combination of the product of projectors on local 
states as in Eq. (j4]). Otherwise, the mixed state is said to be entangled. 

The number of pure separable states K necessary to decompose any separable state 
according to Eq. (j4]) is limited by the Caratheodory theorem as K < [dAcLBY (see 
Refs. |[8 10 1). No better bound is known in general. 



By definition, entangled states cannot be prepared locally by two parties even af- 
ter communicating over a classical channel. To prepare entangled states, the physical 
systems must be brought together to interacj^ Mathematically, a nonlocal unitary 
operatoij^] must necessarily act on the physical system described by T-La <8i to 
produce an entangled state from an initial separable state. 

The question whether a given bipartite state is separable or not turns out to be 
quite complicated. Although the general answer to the separability problem still 
eludes us, there has been significant progress in recent years, and we will review 
some such directions in the following paragraphs. 



2.3 Entanglement criteria 

An operational necessary and sufficient criterion for detecting entanglement still 
does not exist. However, over the years the whole variety of criteria allowing for 
detection of entanglement has been worked out. Below we review some of the most 
important ones, while for others the reader is referred to Ref 1 12|. Note that, even 
if we do not have necessary and sufficient separability criteria, there are numerical 
checks of separability: semidefinite programming was used to show that separability 
can be tested in a finite number of steps, although this number can become too 
large for big systems |13 , 14J. In general — without a restriction on dimensions — 
the separability problem belongs to the NP-hard class of computational complexity 




' Due to entanglement swapping 1 1 1 1, one must suitably enlarge the notion of preparation of en- 
tangled states. So, an entangled state between two particles can be prepared if and only if either 
the two particles (call them A and B) themselves come together to interact at a time in the past, or 
two other particles (call them C and D) do the same, with C having interacted beforehand with A 
and D with B. 

^ A unitary operator on Ha ^Hb is said to be "nonlocal" if it is not of the form Ua Ub , where 
Ua is a unitary operator acting on Ha and Ub acts on "Hb- 
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2.4 Partial Transposition 



Let us start with an easy-to-apply necessary criterion based on the transposition 
map recognized by Choi |16| and then independently formulated directly in the 
separability context by Peres 1 17|. 

Let qab be a state on the product Hilbert space Hab, and T : B{C^) — ^ B{C^) 
a transposition map with respect to the some basis {\i)} in C^, defined through 
T{X) = = J:^.J^^J\M\ for any X = ^^Mj\ from BiC). Let us 
now consider an extended map T ® Ib called hereafter partial transposition, where 
Ib is the identity map acting on the second subsystem. When applied to qab, the 
map T transposes the first subsystem leaving the second one untouched. More 
formally, writing qab as 

cLa — 1 ds — 1 

Q^B=Y. E Qt;m\®\M (5) 

where and {\^)} are bases in Alice and Bob Hilbert spaces, respectively, we 
have 

Na Nb 

{t®ib){qab) ^ Q% = E E CIjX^i ® Ia^xh- (6) 

— 1 1 

Similarly, one may define partial transposition with respect to the Bob's subsys- 
tem (denoted by p^^). Although the partial transposition of qab depends upon the 
choice of the basis in which qab is written, its eigenvalues are basis independent. 
The applicability of the transposition map in the separability problem can be for- 
maUzed by the following statement. 

Theorem 2. 1 17 If a state pab is separable, then > Q and > 0. 



Proof. Since qab is separable, according to definition [2] it has the form Q. Then, 
performing the partial transposition with respect to the first subsystem, we have 



K K 



In the second step we used that ~ (^*)^ for all A. The above shows that 

is a proper (and also separable) density matrix implying that p^'^ > 0. The same 

reasoning leads to the conclusion that > 0, finishing the proof. ■ 

Due to the identity Qj^^ = (g*^^)^, and the fact that global transposition does not 
change eigenvalues, partial transpositions with respect to the A and B subsystems 
are equivalent from the point of view of the separability problem. 

In conclusion, we have a simple criterion {partial transposition criterion) for 
detecting entanglement. More precisely, if the spectrum of one of the partial trans- 
positions of QAB contains at least one negative eigenvalue then qab is entangled. As 
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an example, let us apply the criterion to pure entangled states. If | V'As) is entangled, 
it can be written as ^ with r > 1. Then, the eigenvalues of \'4'ab){'4'ab\'^^ will be 
Af (i = 1, . . . , r) and ±\i\j {i j i,j = 1, . . . , r). So, an entangled I^jab) of 
Schmidt rank r > 1 has partial transposition with r(r — l)/2 negative eigenvalues 
violating the criterion stated in theorem[2] 

The partial transposition criterion allows to detect in a straightforward manner all 
entangled states that have non-positive partial transposition (hereafter called NPT 
states). However, even if this is a large class of states, it turns out that — as pointed 
out in Refs. | T0p8 1 — there exist entangled states with positive partial transposition 



(called PPT states) (cf. Fig. |2]i. Moreover, the set of PPT entangled states does 
not have measure zero | [T9| . It is, therefore, important to have further independent 
criteria that identifies entangled PPT states. Remarkably, PPT entangled states are 
the only known examples of bound entangled states, i.e., states from which one 
cannot distill entanglement by means of LOCC, even if the parties have an access 
to an unlimited number of copies of the state ||8]|T8). The conjecture that there exist 
NPT "bound entangled" states is one of the most challenging open problems in 
quantum information theory pO||2T| . Note also that both separable as well as PPT 
states form convex sets. 

Theorem |2] is a necessary condition of separability in any arbitrary dimension. 
However, for some special cases, the partial transposition criterion is both a neces- 
sary and sufficient condition for separability: 

Theorem 3. ]22) A state qab acting on (E) or (E) is separable if and 
only if 0% > 0- 

We will prove this theorem later. Also, we will see that Theorem|2]is true for a whole 
class of maps (of which the transposition map is only a particular example), which 
also provide a sufficient criterion for separability. Before this, let us discuss the dual 
characterization of separability via entanglement witnesses. 



2.5 Entanglement Witnesses from the Hahn-Banach theorem 

Central to the concept of entanglement witnesses is the corollary from the Hahn- 
Banach theorem (or Hahn-Banach separation theorem), which we will present here 
limited to our needs and without proof (which the reader can find e.g. in Ref. p3j). 

Theorem 4. Let S be a convex compact set in a finite-dimensional Banach space. 
Let p be a point in this space, however, outside the set S (p ^ S). Then there exists 
a hyperplan^that separates pfrom S. 

The statement of the theorem is illustrated in figure [T] In order to apply it to our 
problem let S denote now the set of all separable states acting on T-La®T^b- This is a 
convex compact subset of the Banach space of all the linear operators B{HaE)Hb)- 



' A hyperplane is a linear subspace with dimension one less than the dimension of the space itself. 
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Fig. 1 Schematic picture of the Hahn-Banach theorem. The (unique) unit vector orthonormal to 
the hyperplane can be used to define right and left with respect to the hyperplane by using the sign 
of the scalar product. 

The theorem implies that for any entangled state Qj^b there exists a hyperplane 
separating it from S. 

Let us introduce a coordinate system located within the hyperplane (along with 
an orthogonal vector W chosen so that it points towards S). Then, every state 
gAB can be characterized by its "distance" from the plane, here represented by 
the Hilbert-Schmidt scalar producj^ According to our choice of the coordinate sys- 
tem (see Fig. [TJ, for any such hyperplane W every separable state has a positive 
"distance", while there are some entangled states with a negative "distance". More 
formally, theorem Q implies the following seminal result. 

Theorem 5. ]22) Let qab be some entangled state acting on Has- Then there 
exists a Hermitian operator W G B{T-La ® "Hb) such that tT:{gABW) < and 
tT{c7ABW) > O for all separable (tab & BiHAf^HB)- 

It is then clear that all the operators W representing such separating hyperplanes 
deserve special attention as they are natural candidates for entanglement detec- 
tors. That is, given some Hermitian W, if tT{W qab) < and simultaneously 
ti{WaAB) > for all separable aAB, we know that qab is entangled. One is 
then tempted to introduce the following definition |24| . 

Definition 3. We call the Hermitian operator W an entanglement witness if 
tY{WaAB) > for all separable aAB and there exists an entangled state qab such 
that tr{W qab) < 0. 



Let be some Hilbert space. Then the set B('H) of linear bounded operators acting on "H is also 
a Hilbert space with the Hilbert-Schmidt scalar product {A\B) = tr(At_B) (A, B e 
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Example 1. Let us discuss how to construct entanglement witnesses for all NPT 
states. If qab is NPT then its partial transposition has at least one negative eigen- 
value. Let IV'i) denote the eigenstates of corresponding to its negative eigenval- 
ues Xi < 0. Then the Hermitian operator Wi = has negative mean value 
on qab, i.e., tr(eAB|^/'i)('0ir^) = trCt'AslV'jXV'il) = < 0. Simultaneously, 
using the identity tr{AB'^) = ti{A^B) obeyed by any pair of matrices A and B, 
it is straightforward to verify that ti{Wia'AB) > for all i and separable gab- One 
notices also that any affine combination of Wi and in particular g^^g itself are also 
entanglement witnesses. 

Let us comment shortly on the properties of entanglement witnesses. First, it is 
clear that they have negative eigenvalues, as otherwise their mean value on all entan- 
gled states would be positive. Second, since entanglement witnesses are Hermitian, 
they can be treated as physical observables — which means that separability cri- 
teria based on entanglement witnesses are interesting from the experimental point 
of view. Third, even if conceptually easy, entanglement witnesses depend on states 
in the sense that there exist entangled states that are only detected by different wit- 
nesses. Thus, in principle, the knowledge of all entanglement witnesses is necessary 
to detect all entangled states. 



2.6 Positive maps and the entanglement problem 

Transposition is not the only map that can be used to deal with the separability 
problem. It is rather clear that the statement of theorem [2] remains true if, instead 
of the transposition map, one uses any map that when applied to a positive operator 



gives again a positive operator (a positive map). Remarkably, as shown in Ref. 1 22 1, 
positive maps give not only necessary but also sufficient conditions for separability 
and entanglement detection. Moreover, via the Jamiolkowski-Choi isomorphism, 
theorem |5] can be restated in terms of positive maps. To see this in more detail we 
need to review a bit of terminology. 

We say that a map ^1 : B{Ha) ^ B{Hb) is linear if A{aX + l3Y) = aAiX) + 
PA{Y) for any pair of operators X, Y acting on Ha and complex numbers a, /3. We 
also say that A is Hermiticity-preserving (trace-preserving) if A{X^) — [A(X)]^ 
(tT[A{X)] = tr{X)) for any Hermitian X e B{Ha)- 

Definition 4. A linear map A : B{Ha) B{Hb) is called positive if for all positive 
X e B{Ha) the operator A{X) e B{Hb) is positive. 

As every Hermitian operator can be written as a difference between two positive 
operators, any positive map is also Hermiticity-preserving. On the other hand, a 
positive map does not have to be necessarily trace-preserving. 

It follows immediately from the above definition that positive maps applied to 
density matrices give (usually unnormalized) density matrices. One could then ex- 
pect that positive maps are sufficient to describe all quantum operations (as for 
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instance measurements). This, however, is not enough, as it may happen that the 
considered system is only part of a larger one and we must require that any quan- 
tum operation on our system leaves the global system in a valid physical state. This 
requirement leads us to the notion of completely positive maps: 

Definition 5. Let A : BCHa) ~^ B{'Hb) be a positive map and let Id : Md{C) 
AId{C) denote an identity map. Then, we say that A is completely positive if for all 
d the extended map Id<E} Ais positive. 

Let us illustrate the above definitions with some examples. 

Example 2. (Hamiltonian evolution of a quantum state) Let Ha = Hb = H and 
let Au : B{H) B{H) be defined as Au{X) = UXW for any X e B{H), 
with U being some unitary operation acting on Ti. Since unitary operations do not 
change eigenvalues when applied to X, it is clear that Au is positive for any such U. 
Furthermore, Au is completely positive: an application of the extended map Id^Ajj 
toX e BiH-^H) gives {Id(8)Au)iX) = {ld(g>U)X{ld'S>U)\ where l<j denotes 
identity acting on T-L. Therefore, the extended unitary U — Id <E) U is also unitary. 
Thus, if X > 0, then UXW > 0. The commonly known example of Ajj is the 
unitary evolution of a quantum state g{t) = U (t) g{0)[U (t)]^ = yl(7(f)(gi(0)). 

Example 3. (Transposition map) The second example of a linear map is the already 
considered transposition map T. It is easy to check that T is Hermiticity and trace- 
preserving. However, the previously discussed example of partially transposed pure 
entangled states shows that it cannot be completely positive. 

To complete the characterization of positive and completely positive maps let 
us just mention the Choi-Kraus-Stinespring representation. Recall first that any 
linear Hermiticity-preserving (and so positive) map A : B{C^) B{C^) can be 
represented as f25l : 

k 

-^?7,y,Xt//, (8) 
1=1 

where k < <P, rji £ M, and V^; : — > C'' are orthogonal in the Hilbert-Schmidt 
scalar product tr(y/Vj) — Sij. In this representation, completely positive maps are 
those (and only those) that have rji > for all i. As a result, by replacing Wi = 
^Jffi Vi (which preserves the orthogonality of Wi), we arrive at the aforementioned 
form for completely positive matrices p6[|28| . 

Tlieorem 6. A linear map A : B{C^) — > B{C^) is completely positive iff admits the 
Choi-Kraus-Stinespring fon7i 

k 

A{X)^J2^^^^^^ (9) 

i=l 

where k < d^ and ; — > C^, called usually Kraus operators, are orthogonal 
in the Hilbert-Schmidt scalar product. 
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Finally, let us recall the so-called Choi-Jamiolkowski isomorphism p6l[29) : ev- 
ery Unear operator X acting on (X) can be represented as X = [I ® A) (P!f ^ ) 
with some linear map A : B{C'^) — > B{C^). With this isomorphism, entangle- 
ment witnesses correspond to positive maps. Notice also that the dual form of this 
isomorphism reads A{X) = tTB[W{lA <8i X'^)])- 

Equipped with new definitions and theorems, we can now continue with the rela- 
tionship between positive maps and the separability problem. It should be clear by 
now that theorem |2] is just a special case of a more general necessary condition for 
separability: if qab acting on Ha ® Hb is separable, then (/ (g) A){gAB) is positive 
for any positive map A. In a seminal paper in 1996 |22|, Horodecki and Horodecki 
showed that positive maps also give a sufficient condition for separability. More 
precisely, they proved the following: 

Theorem?. |22j A state pab G B{'C'^^ ® C^^) is separable if and only if the 
condition 

{I (E> A)ipAB) > 0. (10) 
holds for for all positive maps A : 6(C''«) S(C'''*). 

Proof. The "only if part goes along exactly the same lines as proof of theorem |2j 
where instead of the transposition map we put A. 

On the other hand, the "if part is much more involved. Assuming that qab is 
entangled, we first show that there exists a positive map A : B{C^^) —?' B{C^^) 
such that (/ (X) A){gAB) ^ 0. For this we can use theorem|5] which says that for 
any entangled qab there always exists entanglement witness W detecting it, i.e., 
ti{WgAB ) < 0. Denoting by L : B{C^^ ) — ^ S(C'^^ ) a positive map corresponding 
to the witness W via the the Choi-Jamiolkowski isomorphism, i.e., W — {I ^ 
L){p'f^^), we can rewrite this condition as 

ti[{I(g>L){Pi^^^)gAB]<0. (11) 

As L is positive it can be represented as in Eq. ([8]), and hence the above may be 
rewritten as Tr[pf^\l®L^){gAB)] with : 6(C'^«) ^ BiC^^) called the dual 
map of L. One immediately checks that dual maps of positive maps are positive. 
This actually finishes the proof since we showed that there exists a positive map 
A = L'' such that (/ (g) A){gAB) ^ 0. ■ 

In conclusion, we have two equivalent characterizations of separability in bipartite 
systems, in terms of either entanglement witnesses or positive maps. However, on 
the level of a particular entanglement witness and the corresponding map, both char- 
acterizations are no longer equivalent. This is because usually maps are stronger in 
detection than entanglement witnesses (see Ref. |30|). A good example comes from 
the two qubit case. On one hand, theorem[3]tell us that the transposition map detects 
all the two-qubit entangled states. On the other hand, it is clear that the correspond- 
ing witness, the so-called swap operator (see Ref. l9|) V = pj^^^ does not detect 
all entangled states — as for instance tr(P| 'V) > 0. 
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Let us also notice that an analogous theorem was proven in Ref. pO| , which gave 
a characterization of the set of the fully separable multipartite states 

qa,...Am =^PiQ'm® S% (12) 

i 

in terms of multipartite entanglement witnesses. Here, however, instead of positive 
maps one deals with maps which are positive on products of positive operators. 



2.7 Positive maps and entanglement witnesses: further 
characterization and examples 

We discuss here the relationship between positive maps (or the equivalent entangle- 
ment witnesses) and the separability problem. 

Deflnition 6. Let A : B{Ha) ~> ^(^b) be a positive map. We call it decomposable 
if it admits the forrrj^ = Af^ + A^^ o T, where Af^ {i = 1,2) are some 
completely positive maps. Otherwise A is called indecomposable. 

It follows from this definition that decomposable maps are useless for detection 
of PPT entangled states. To see this explicitly, assume that qab is PPT entangled. 
Then it holds that (/ ^ A){gAB) = (/ A^P){gAB) + {I <E> ^?^)((?Il3) = (/ «) 
■^i^){0ab) + (/ (X" A2^){£ab), where qab = Sab some quantum state. Since 
Af^ are completely positive, both terms are positive and thus (/ ® A){qab) > 
for any decomposable A and PPT entangled qab ■ 

The simplest example of a decomposable map is the transposition map, with both 
Af^ [i — 1,2) being just the identity map. It is then clear that, from the point of 
view of entanglement detection, the transposition map is also the most powerful 
example of a decomposable map. Furthermore, as shown by Woronowicz pT[ , all 
positive maps from B{€?) and B{'C^) to B{€?) are decomposable. Therefore, the 
partial transposition criterion is necessary and sufficient in two-qubit and qubit- 
qutrit systems as stated in theorem|3] 

Using the Jamiolkowski-Choi isomorphism we can check the form of entangle- 
ment witnesses corresponding to the decomposable positive maps. One immediately 
sees that they can be written as = P + , with P and Q being some positive 
operators. Following the nomenclature of positive maps, such witnesses are called 
decomposable. 

It is then clear that PPT entangled states can only be detected by indecomposable 
maps, or, equivalently indecomposable entanglement witnesses (cf. Fig. |2]). Still, 
however, there is no criterion that allows to judge unambiguously if a given PPT 
state is entangled. 



^ By yli o A2 we denote the composition of two maps Ai (i = 1, 2), i.e., a map that acts on a given 
operator X as yli o A2iX) = Ai(A2(X)). 
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Fig. 2 Schematic view of the Hilbert-space with two states pi and p2 and two witnesses EWl 
and EW2. EWl is a decomposable EW, and it detects only NPT states like pi. EW2 is an 
indecomposable EW, and it detects also some PPT states like p2 . Note that none of the witnesses 
detect all entangled states. 



To support the above discussion, we give particular examples of positive maps 
and corresponding entanglement witnesses. 

Example 4. Let Ar : ;B(C'^) — > ,B(C'^) be the so-called reduction map map defined 
through Ar{X) = tr(X)lrf - X for any X £ S(C''). It was introduced in Ref. |[32) 
and considered firstly in the entanglement context in Refs. 1 33 34 1 . One immediately 
finds that Ar is positive but not completely positive, as it detects entanglement of 
pI^\ Moreover, Ar = A'^^oT, where yl^P is a completely positive map with Kraus 
operators (cf. theorem|6]) given by Vij ~ — {i < j, i, i = 0, . . . , d — 1), 
meaning that the reduction map is decomposable. 

Example 5. Let A^^^ : BIC^) — > ,B(C'^) be the so-called extended reduction map 
|[35,36|definedbyyl^^t(X) = tr(X)ld-X- ^7X^^7^ where [/ obeys U'^ = -U 
wdlPU < Id At is obviously positive but not completely positive. However, unlike 
the reduction map, this one is indecomposable as examples of PPT entangled states 
detected by yl^j can be found |35 , 36 1. 

Let us summarize our considerations with the following two theorems. First, us- 
ing the definitions of decomposable and indecomposable entanglement witnesses, 
we can restate the consequences of the Hahn-Banach theorem in several ways: 



Theorem 8. 1 16 22 37 -39) The following statements hold. 



1. A state pab is entangled iff there exists an entanglement witness W such that 

tl{WpAB) < 0. 

2. A state pab is PPT entangled iff there exists an indecomposable entanglement 
witness W such that tT{W pab) < 0- 

3. A state uab separable iff tT{WaAB) ^ 0/or all entanglement witnesses. 
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Notice that the JamioBcowski-Choi isomorphism between positive maps and entan- 
glement witnesses allows to rewrite immediately the above theorem in terms of posi- 
tive maps. From a theoretical point of view, the theorem is quite powerful. However, 
it does not give any insight on how to construct for a given state p, the appropriate 
witness operator. 

Second, the relations between maps and witnesses can be collected as follows. 

Theorem 9. p2||29|[37] - [39[ Let W be a Hennitian operator and A\y map defined 
as Aw{X) = trB[W{lA ® X'^)]. Then the following statements hold. 

1. W >Q iff Aw is a completely positive map. 

2. W is an entanglement witness iff Aw is a positive map. 

3. W is a decomposable entanglement witness iff Aw is decomposable map. 



2.8 Entanglement measures 



The criteria discussed above allow to check if a given state qab is entangled. How- 
ever, in general they do not tell us directly how much qab is entangled. In what 
foUows we discuss several methods to quantify entanglement of bipartite states. 
This quantification is necessary, at least partly because entanglement is viewed as a 
resource in quantum information theory. There are several complementary ways to 
quantify entanglement (see Refs. |[8] |4QH49| and references therein). We will present 
here three possible ways to do so. 

Let us just say few words about the definition of entanglement measure^ The 
main ingredient in this definition is the monotonicity under LOCC operations. More 
precisely, if A denotes some LOCC operation, and E our candidate for the entan- 
glement measure, E has to satisfy 

E{A{q)) < E{e) (13) 



or 



where Qi are states resulting from the LOCC operation A appearing with proba- 
bilities Pi (as in the case of e.g. projective measurements). Both requirements fol- 
low from the very intuitive condition saying that entanglement should not increase 
under local operations and classical communication. It follows also that if E is 



convex, then the condition ( 14 1 implies ( 13 1, but not vice versa — therefore ( 14 1 



gives a stronger condition for the monotonicity. For instance, the three examples of 
measures presented below satisfy this condition. Finally, notice that from the mono- 
tonicity under LOCC operations one also concludes that E is invariant under unitary 
operations, and gives a constant value on separable states (see e.g. Ref. |[8|). 

* For a more detailed axiomatic description, and other properties of entanglement measures, the 
reader is encourae ed to consult, e.g., Refs. (8]|48[49) . 
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2.8.1 Entanglement of formation 

Consider a bipartite pure state IV'ab) & C'*'* C''^ shared between Alice and 
Bob. As shown by Bennett et al. pO) , given nE{\ijjAB)) copies of the maximally 
entangled state, Alice and Bob can by LOCC transform them into n copies of | ipAs) , 
if n is large. Here 

E{\^ab)) - SiQA) - Sige) (15) 

with QA and qb being the local density matrices of \iPab) and S{p) stands for 
the von Neumann entropy of p given by S{p) — -~trp\og2 p. It clearly follows 
from theorem [T] that E is zero iff \tpAB) is separable, while its maximal value 
log2 mm{dA,dB} is attained for the maximally entangled states [\\. 

(2) — 

For the two-qubit maximally entangled state ), the function E gives one: 
an amount of entanglement also called ebit. With this terminology, one can say that 
\tpAB) has E{\iPab)) ebits. Since E{iPab) is the number of singlets required to 
prepare a copy of the state \iPab), it is called entanglement of formation of \4'ab)- 
We are therefore using the amount of entanglement of the singlet state as our unit of 
entanglement. 



Following Ref. |40|, let us now extend the definition of entanglement of forma- 
tion to all bipartite states. By definition, any mixed state is a convex combination of 
pure states, i.e., g — pi\%l]i){'4)i\, where probabilities pi and pure states (not nec- 
essarily orthogonal) constitute what is called an ensemble. A particular example 
of such an ensemble is the eigendecomposition of g. Thus, it could be tempting to 
define the entanglement of formation of g as an averaged cost of producing pure 
states from the ensemble, i.e., X^i Pi^dV'i))- One knows, however, that there exist 
an infinite number of ensembles realizing any given g. A natural solution is then to 
minimize the above function over all such ensembles — with which we arrive at the 



definition of entanglement of formation for mixed states pO) : 

E{gAB)^ mill Y^p,E{\i{;,)), (16) 

with the minimum taken over all ensembles {p^, such that X^iP^IV'iXV'il — 
Qab- 

In general, the above minimization makes the calculation of entanglement of 



formation extremely difficult. Nevertheless, it was determined for two-qubits 1 5 1 
[S?!, or states having some symmetries, as the so-called isotropic p3) and Werner 
|54j| states. In the first case it amounts to 

p , . ^i ^ + ^^-C\gAB) \ 

EpigAB) = H , (17) 



where H{x) — —x logj x — {1 — x) log2(l — x) is the binary entropy function. The 
function C is given by 



C{gAB) = max{0, Ai - A2 - A3 - A4} . 



(18) 
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withAi, . . . , A4the eigenvalues of the Hermitian matrix [{qabY^'^ Qab^QabY^'^Y^'^ 
in decreasing order, and qab = (Jy <E) o'y9*AB'^y ® '^y Note that the complex con- 
jugation over g is taken in the cr^ eigenbasis, and ay denotes the well-known Pauli 
matrij(|^ The function C, called concurrence, can also be used to quantify entangle- 
ment of more general quantum states. Although Eq. (18 1 gives the explicit form of 
concurrence only for two-qubit states, it can also be defined for arbitrary bipartite 
states — as we shall discuss in the following section. 



2.8.2 Concurrence 



For any \4'ab) & C^'^'S^C^" we define concurrence as CQiPab)) — ^2(1 — trp^) 
where is one of the subsystems of \4>ab) (note that the value of C does not 
depend on the choice of subsystems) | |55) . In the case dA = ds = d, one sees that 
its value for pure states ranges from for separable states to y^2{l — 1/d) for the 
maximally entangled state. 

The extension to mixed states goes in exactly the same way as in the case of 
entanglement of formation, 

C(e^B)= min (19) 

{P^^i',)} 

I 

where again the minimization is taken over all the ensembles that realize gAB- For 
the same reason, as in the case of EOF, concurrence is calculated only in few in- 
stances like two-qubit states |51||52| and isotropic states p6) . 

Seemingly, the only difference between E and C lies in the function taken to de- 
fine both measures for pure states. However, the way concurrence is defined enables 
one to determine it experimentally for pure states 1 57 58), provided that two copies 
of the state are available simultaneously. 



2.8.3 Negativity and logarithmic negativity 

Based on the previous examples of entanglement measures, one may get the impres- 
sion that all of them are difficult to determine. Even if this is true in general, there 
are entanglement measures that can be calculated for arbitrary states. The examples 
we present here are negativity and logarithmic negativity. The first one is defined 
as KT91I591: 

N{gAB) = \{Wg'AB\\-^)- (20) 

The calculation of N even for mixed states reduces to determination of eigenvalues 
of g^%, and amounts to the sum of the absolute values of negative eigenvalues of 
g^B- This measure has a disadvantage: partial transposition does not detect PPT 



^ In the standard basis ay is given by ay = 



-i|0>(l| + i|l>(0|. 
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entangled states; therefore N is zero not only for separable states but also for all 
PPT states. 

The logarithmic negativity is defined as p9) : 

ENiQAB) = l0g2 WqabW = l0g2[2iV(eAB) + 1] ■ (21) 

It was shown in Ref. | [60) that it satisfies condition ( [T4| . Moreover, logarithmic 
negativity is additive, i.e., E{qab ® <^ab) — E{qab) + E{aAB) for any pair of 
density matrices qab and gab, which is a desirable feature. However, this comes 
at a cost: Ejq is not convex |60|. Furthermore, for the same reason as negativity 
it cannot be used to quantity entanglement of PPT entangled states. Finally, let us 
notice that these measures range from zero for separable states, to [d — l)/2 for 
negativity and log2 d for logarithmic negativity. 



3 Area laws 



Area laws play a very important role in many areas of physics, since generically 
relevant states of physical systems described by local Hamiltonians (both quantum 
and classical) fulfill them. This goes back to the seminal work on the free Klein- 
Gordon field | 6T||62) , where it was suggested that the area law of geometric entropy 



might be related to the physics of black holes, and in particular the Bekenstein- 
Hawking entropy that is proportional to the area of the black hole surface [63 - 65| . 
The related holographic principle | [66| says that information about a region of space 
can be represented by a theory which lives on a boundary of that region. In recent 
years there has been a wealth of studies of area laws, and there are excellent reviews 
| [67) and special issues |68| about the subject. As pointed out by the authors of 
Ref. | [67) , the interest in area laws is particularly motivated by the four following 
issues: 

• The holographic principle and the entropy of black holes, 

• Quantum correlations in many-body systems, 

• Computational complexity of quantum many-body systems, 

• Topological entanglement entropy as an indicator of topological order in certain 
many-body systems 



3.1 Mean entanglement of bipartite states 

Before we turn to the area laws for physically relevant states let us first consider a 
generic pure state in the Hilbert space in C™ (g) C" (to < n). Such a generic state 
(normalized, i.e. unit vector) has the form 
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l'^)-EE"^^l*)l^')' (22) 

i=l 3 = 1 

where the complex numbers aij may be regarded as random variables distributed 
uniformly on a hypersphere, i.e. distributed according to the probability density 

(m n \ 

with the only constraint being the normalization. As we shall see, such a generic 
state fulfills on average a "volume" rather than an area law. To this aim we in- 
troduce a somewhat more rigorous description, and we prove that on average, the 
entropy of one of subsystems of bipartite pure states in C™ (E) C" (m < n) is almost 
maximal for sufficiently large n. In other words, typical pure states in C" (S) C" 
are almost maximally entangled. This "typical behavior" of pure states happens to 
be completely atypical for ground states of local Hamiltonians with an energy gap 
between ground and first excited eigenstates. 



Rigorously speaking, the average with respect to the distribution ( 23 1 should 
be taken with respect to the unitarily invariant measure on the projective space 
(j-^p?riri-i jj. ^ unique measure generated by the Haar measure on the unitary group 
by applying the unitary group on an arbitrarily chosen pure state. One can show then 
that the eigenvalues of the first subsystem of a randomly generated pure state \iPab) 
are distributed according to the following probability distribution 1 69 ■ 71] (see also 
Ref. my. 



P,nA^i, . . . , A™) = C„,,„(5(^ A, - 1) n Ar™ [](A, - A,)2, (24) 

i i i<j 

where the delta function is responsible for the normalization, and the normalization 
constant reads (see e.g. Ref. | |72j ) 

C„ „ = , (25) 

with 7^ being the Euler gamma functiorj^ 

Theorem 10. Lef \iPab) be a bipartite pure state from C" ® C" (m < n) 
drawn at random according to the Haar measure on the unitary group and qa = 

* In general the gamma function is defined through 

r{z) = I e-^e-^dt (z e C). (26) 





For z being positive integers z = n the gamma function is related to the factorial function via 
r(n) = (n - 1)!. 
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tTB\'ipAB){'4'AB \ be its subsystem acting on C™. Then, 

TTl 

{S{gA))~\ogm-—. (27) 
in 

Proof. Let us give here just an intuitive proof without detailed mathematical discus- 
sion (which can be found e.g. in Refs. f69^-'75l). 
Our aim is to estimate the following quantity 



{S{qa)) 



J log A,^ P(Ai, . . . , A™) dAi . . . dAi, (28) 



where the probability distribution P(Ai, . . . , Am) is given by Eq. (24 1. We can al- 
ways write the eigenvalues Ai as A; = — + Si, where Si £ M and - Si — 0. This 
allows us to expand the logarithm into the Taylor series in the neighborhood of 1/m 
as 

log ( - + 5, ) = - log m + ^ ^ (mS^ f, (29) 



k=l 



which after application to Eq. ( 28 i gives the following expression for the mean 
entropy 

(30) 

Let us now notice that tTg\ = A| = ^i{Si + = Sf + 1/m, and 

therefore J2i — ^"^Qa — This, after substitution in the above expression, 
together with the fact that for sufficiently large n we can omit terms with higher 
powers of Si (cf. |[69)), leads us to 

{S{QA))^\ogm-'^(ire\-^^. (31) 

One knows that irg^ denotes the purity of qa. Its average was calculated by Lubkin 
and reads 

/+ 2\ m + n 
^ ' mn + 1 



Substitution in Eq. ( 3 1 1 leads to the desired results, completing the proof. ■ 

Two remarks should be made before discussing area laws. First, it should be 
pointed out that it is possible to get analytically the exact value of (5). There is a 
series of papers |j73}|75l presenting different approaches leading to 



771 — 1 

{S{q)) ^ ^{mn + 1) - ^{n +1) (33) 

Zn 
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with if' the bigamma functiorj^ Using now the fact that ^{z + 1) = + 1/z, 
and the asymptotic properties of bigamma function, ^'(z) ~ log z, we get (27 1. 

Second, notice that the exact resuh of Lubkin ( [32| can be estimated by relaxing 
the normalization constraint in the distribution (|23|, and replacing it by a product of 
independent Gaussian distributions, P{a) — Yli jinm/ir) exp[—nm\aij\'^], with 
(aij) — 0, and p) = 1/nm. The latter distribution, according to the central 
limit theorem, tends for nm — > oo to a Gaussian distribution for l^^ij P 

centered at 1, with width ~ 1/y/nm. One obtains then straightforwardly (trgA) = 
1, and after a little more tedious calculation {trg\) = {n + m)/nm, which agrees 
asymptotically with the Lubkin result for nm ^ 1. 



3.2 Area laws in a nutshell 

In what follows we shall be concerned with lattices L in D spatial dimensions, 
L C Z^. At each site we have a d-dimensional physical quantum system (one can, 
however, consider also classical lattices, with a rf-dimensional classical spin at each 
site with the configuration space = {0, . . . ,d — 1}) at each sit^^ The distance 
between two sites x and y of the lattice is defined as 

'D{x,y)^ max \x,-y.^. (35) 

l<i<D 

Accordingly, we define the distance between two disjoint regions X and F of L as 
the minimal distance between all pairs of sites {x,y}, where x £ X and y ^ Y; 
i.e., 'D{X,Y) — miiijjgx minjjgy y). If R is some region of L, we define 
its boundary dR as the set of sites belonging to R whose distance to L \ R (the 
complement of R) is one. Formally, dR — {x ^ R\'D{x, L \ R) = 1}. Finally, by 
\R\ we denote number of sites (or volume) in the region R (see Figure [3]). 

Now, we can add some physics to our lattice by assuming that interactions be- 
tween the sites of L are governed by some hamiltonian H. We can divide the lattice 
L into two parts, the region R and its complement L\R. Roughly speaking, we aim 
to understand how the entropy of the subsystem R scales with its size. In particular, 
we are interested in the entropy of the state grt reduced from a ground state or a 
thermal state of the Hamiltonian H. We say that the entropy satisfies an area law if 



' The bigamma function is defined as 'I'{z) = F' (z)/ r(z) and for natural z = n it takes the form 



n-l 



1'{n) = -7 + y" - (34) 

z — ^ r) 



fc = l 



with 7 being the Euler constant, of which exact value is not necessary for our consideration as it 
vanishes in Eq. l |33[ l. 

10 pqj. results concerning other kind of systems one can consult Ref . |67| . 
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Fig. 3 Schematic representation of a lattice system L, an arbitrary region R (denoted in light grey 
background), and its boundary dR (denoted in dark grey background). 



it scales at most as the boundary area'M i.e., 



3.2.1 One-dimensional systems 

Let us start with the simplest case of one-dimensional lattices, L = {1, . . . , iV}. Let 
R be a subset of L consisting of n contiguous spins starting from the first site, i.e., 
R — {\, . . . ,n} with n < N. In this case the boundary dR of the region R contains 
one spin for open boundary conditions, and two for periodic ones. Therefore, in this 
case the area law is extremely simple: 



The case of 13 = 1 seems to be quite well understood. In general, all local 
gapped systems (away from criticality) satisfy the above law, and there might be 
a logarithmic divergence of entanglement entropy when the system is critical. To 
be more precise, let us recall the theorem of Hastings leading to the first of the 
above statements, followed by examples of critical systems showing a logarithmic 
divergence of the entropy with the size of R. 

Consider the nearest-neighbor interaction Hamiltonian 



Sign) = 0{\dR\). 



(36) 



(37) 




(38) 



" Let us shortly recall that the notation f{x) = 0{g{x)) means that there exist a positive constant 
c and xq > such that for any x > a;o it holds that f{x) < cg{x). 
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where each has a nontrivial support only on the sites i and i + 1. We assume 



also that the operator norm of all the terms in Eq. (38 i are upper bounded by some 
positive constant J, i.e., < J for all i (i.e., we assume that the interaction 

strength between ith site and its nearest-neighbor is not greater that some constant). 



Under these assumptions, Hastings proved ]76) the following: 

Theorem 11. Let L be a one-dimensional lattice with N d-dimensional sites, and 



let H be a local Hamiltonian as in Eq. ( 38 ). Assuming that H has a unique ground 



state separated from the first excited states by the energy gap AE > 0, the entropy 
of any region R satisfies 

<6coC2^«'°s'' lege log d (39) 

with Co denoting some constant of order unity and ^ — mm{2v /AE,^c}- Here, 
V denotes the sound velocity and is of order J, while is a length scale of order 
unity. 

Let us remark that both constants appearing in the above theorem come from 
the Lieb-Robinson bound 1 ,77] (see also Ref. [78J for a recent simple proof of this 
bound). 

This theorem tells us that when the one-dimensional system with the local in- 



teraction defined by Eq. (38 i is away from the criticality {AE > 0), the entropy 
of R is bounded by some constant independent of — even if this bound does 
not have to be tight. Of course, we can naturally ask if there exist gapped systems 



with long-range interaction violating (37 1. This was answered in the affirmative in 



Ref. 1 79 80), which gave examples of one-dimensional models with long-range in- 
teractions, nonzero energy gap, and scaling of entropy diverging logaritmically with 
n. 

The second question one could pose is about the behavior of the entropy when 
the gap AE goes to zero and the system becomes critical. Numerous analytical and 
numerical results show that usually one observes a logarithmic divergence of S{gii) 
with the size of the region R. Here we recall only the results obtained for the so- 
called XY model in a transverse magnetic field (for the remaining ones we refer the 
reader to recent reviews ||67][8T), and to the special issue of J. Phys. A devoted to 



this subject ||68|). 

The Hamiltonian for the XY model reads 

ieL ^ ^ ieL 

with < 7 < 1 the anisotropy parameter, and h the magnetic field. In the case 
of vanishing anisotropy (7 = 0), we have the isotropic XY model called shortly 
XX model, while for 7=1 one recovers the well-known Ising Hamiltonian in a 
transverse field. The Hamiltonian Hxy is critical when either 7 = and \h\ < 1 
(the critical XX model) or for \h\ = 1. 

It was shown in a series of papers 1 82 ■ 85] that for the critical X Y model (that is 



when 7 7^ and \h\ = 1) the entropy of the region i? = {1, . . . , n} scales as 
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Sign) = llog^n + 0(1), (41) 

while for the critical XX model, the constant multiplying the logarithms becomes 
one-third. Then, in the case of the critical Ising model (7 = 1), it can be shown that 



the entropy scales at least logaritmicalluH i.e., S{gR) = i7(log2 n) |67 86). 



Concluding, let us mention that there is an extensive literature on the logarithmic 
scaling of the block entropy using conformal field theory methods (see Ref. |j87) 
for a very good overview of these results). Quite generally, the block entropy at 
criticality scales as 

S{gR)^l\og,(^-^^+0{l), (42) 
or, more in general for the Renyi entropjj^ 

Saign) = (c/6) (1 + l/a) log2(|i?|/a) + 0(1), (44) 

where c is called the central charge of the underlying conformal field theory, and a 
is the cutoff parameter (the lattice constant for lattice systems). 



3.2.2 Higher-dimensional systems 

The situation is much more complex in higher spatial dimensions [D > 1). The 



boundary dR of the general area law, Eq. ( 36 1, is no longer a simple one or two 



element set and can have a rather complicated structure. Even if there are no general 



rules discovered so far, it is rather believed that ( 36 1 holds for ground states of lo- 
cal gapped Hamiltonians. This intuition is supported by results showing that for 
quadratic quasifree fermionic and bosonic lattices the area law ( [36| holds |67|. Fur- 
thermore, for critical fermions the entropy of a cubic region R ~ {1, . . . , n}^ is 
bounded as jin^~^ log2 n < S'(pfl) < 72»T.^^^(log2 n)^ with ji [i = 1, 2) denot- 
ing some constants p8]^90) . Let us notice that the proof of this relies on the fact that 
logarithmic negativity (see Eq. ( [2T| i) upper bounds the von Neumann entropy, i.e., 
for any pure bipartite state iV'As), the inequality S{gA{B)) ^ £'Ar(|'0As)) holds. 
This in turn is a consequence of monotonicity of the Renyi entropy with respect 
to the order a, i.e., Sa < S^' for a > a' . This is one of the instances where insights 
from quantum information help to deal with problems in many-body physics. 



The notation f(x) = Q(g{x)) means that there exist c > and xq > such that f{x) > cg{x) 
for all x> xq. 

Recall that the quantum Renyi entropy is defined as 

Sc. = -^log2[Tr(e")] (43) 
I — a 

where a e [0, 00]. For « = one has So{q) = log2 rank(e) and Soo = — log2 Amax with Amax 
being the maximal eigenvalue of q. 
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Interestingly, very recently Masanes (78) showed that the ground state (and also 
low-energy eigenstates) entropy of a region R (even a disjoint one) always scales 
at most as the size of the boundary of R with some correction proportional to 
(log \R\)^ — as long as the Hamiltonian H is of the local form 

H = Y,H,, (45) 

ieL 

where each Hi has nontrivial support only on the nearest-neighbors of the ith site, 
and satisfies as previously \\Hi\\ < J for some J > 0. Thus, the behavior of entropy 
which is considered to be a violation of the area law ([36| can in fact be treated as an 
area law itself. This is because in this cas^^[|9i?|(log — > Ofor|i?| — > oo 
with some A; > 0, meaning that still this behavior of entropy is very different from 
the typical behavior following from theorem [lo] That is, putting m = d'^' and 
n = d)^"^^^ with \L\ ^ \R\ one has that S{gji)/\R.\ is arbitrarily close to logd for 
large \R\. 

Let i?i and i?2 be two disjoint regions of the lattice such that |i?i| < |i?2|, and 
let I denote the distance between these regions. Let us call F a function that bounds 
from above the correlations between two operators X and F(||X||,||F|| < 1) acting 
respectively on i?i and i?2, i.e., C{X,Y) = \{XY) - {X){Y)\ < r{l, \Ri\). The 
first assumption leading to the results of Ref |78| is that if the mean values in C are 
taken in the ground state of H, F is given by 

= ci(Z-Clog|i?i|)"^ (46) 

with some constants ci, ^, and fi > D. Notice that this function decays polynomially 
in I, meaning that this first assumption is weaker than the property of exponential 
decay observed in Ref (91) for gapped Hamiltonians. 

Let now Hj^ denote a part of the global Hamiltonian H which acts only on sites 
in some region R. It has its own eigenvalues and eigenstates, denoted by e„ and 
\ipn) respectively, with eo denoting the lowest eigenvalue. The second assumption 
made in Ref. [78J is that there exist constants C2, t, 7, and 77 such that for any region 
R and energy e = 2J3^\dR\ + eo + 40v, the number of eigenenergies of Hji lower 
than e is upper bounded as 

^Rie) < C2(T|i?|)'^(^-^'')+''l^^l. (47) 

Now, we are in position to formulate the main result of Ref. (78) . 

Theorem 12. Let R be some arbitrary ( even disjoint) region of L. Then, provided 
the assumptions \46\ and \47'^ hold, the entropy of the reduced density matrix gn of 
the ground state of H satisfies 



It should be noticed that one can have much stronger condition for such scahng of entropy. To 
see this expHcitly, say that i? is a cubic region _R = {1, . . . , n}^ meaning that \dR\ = n^~^ and 
|_R| = . Then since Um„_>oo[(logrt)/n'^] = for any (even arbitrarily small) e > 0, one easily 
checks that S{eR)/\dR\^+^ for \dR\ 00. 
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Sign) < C|5i?|(10Clog \R\ f + 0(|ai?|(log \R\ f'^), (48) 

where C collects the constants D, ^, 7, J, 77, and d. If i? is a cubic region, the above 
statement simplifies, giving S{gu) < C\dR\ log \R\ + 0{\dR\) with C being some 
constant. 

Leaving out the first assumption, however, at the cost of extending the second 
assumption to all energies e (not only the ones bounded by 2J3^\dR\ + eo + 401;), 
leads to the following simple area law: 

Theorem 13. Let R be an arbitrary region of the lattice L. Assuming that the above 
number of eigenvalues QR{e) satisfies condition MTvfor all e, then 



S{QR)<C\dR\\og\R\+0{\dR\). (49) 

Proof. Let j^^) and denote the eigenvectors and ordered eigenvalues (cq < ci < 
■ • ■ < e„ < . . .) of Hr. Then, it is clear that the ground state |!fo) of H can be 
written as jtf'o) = Tlii j where the vectors \tfj) constitute some basis 

in the Hilbert space corresponding to the region L\ R. One may always denote 
v^l^i) = I^j "iilV'j)' ™d then 

I'^o) =^Va^|^^)|^.), (50) 

i 

where /k, = = l/X), lo^uP ^ ^^'^ '^hey add up to unity. The vec- 



tors \Lpi) in general do not have to be orthogonal, therefore Eq. (50 1 should not the 
confused with the Schmidt decomposition of |!f'o) ■ Nevertheless, one may show that 
tracing out the L \ R subsystem the entropy of the density matrix acting on R is 
upper bounded as (see Ref. ||50J) 

S{qb) <-^^I.^\0g^i^■ (51) 

i 

We now aim to maximize the right-hand side of the above equation under the fol- 
lowing conditions imposed on First, the locality of our Hamiltonian means that 
{Hr) < gq + J'i^\dR\, implying that the probabilities ^1 obey 

V/^^e, < J3^|9i?|, (52) 



with Bi = Ei — cq. Second, the modified version of the second assumption allows 
to infer that for any eigenvalues the inequality i < C2{T\R\y^'^^''^^^^^ holds. 



Substitution of the above in Eq. ( 52 1 gives 



log^ ^ C-ISi?! log \R\ + 0{\dR\) 

i 



(53) 



26 



R. Augusiak, F. M. Cucchietti, and M. Lewenstein 



where C contains the constants t],j,J, and D. Eventually, following the standard 
convex optimization method (see e.g. Ref. | [92) ) with two constraints (normalization 
and the inequality (53 i) one gets (49 1. ■ 



3.2.3 Are laws for mutual information - classical and quantum Gibbs states 

So far, we considered area laws only for ground states of local Hamiltonians. In 
addition, it would be very interesting to ask similar questions for nonzero tempera- 
tures. Here, however, one cannot rely on the entropy of a ubsystem, as in the case of 
mixed states it looses its meaning. A very good quantity measuring the total amount 
of correlation in bipartite quantum systems is the quantum mutual information | j93j 
defined as 

I{A:B) = S{qa) + S{qb)-S{qab), (54) 

where qab is some bipartite state with its subsystems qa{b)- It should be noticed 
that for pure states the mutual information reduces to twice the amount of entangle- 
ment of the state. 

Recently, it was proven that thermal states gp = e~^^ /\.Y[e~^^\ with local 
Hamiltonians H obey an area law for mutual information. Interestingly, a similar 
conclusion was drawn for classical lattices, in which at each site we have a classical 
spin with the configuration space Z^;, and instead of density matrices one deals with 
probability distributions. In the following we review these two results, starting from 
the classical case. 

To quantify correlations in classical systems, we use the classical mutual infor- 
mation, defined as in Eq. (j54| with the von Neumann entropy substituted by the 
Shannon entropy H{X) = — '}2,,j.p{x) \og2P{x), where p stands for a probability 
distribution characterizing random variable X. More precisely, let A and B = S\A 
denote two subsystems of some classical physical system S. Then, let p{xa) and 
p{xb) be the marginals of the joint probability distribution p{xab) describing S 
(Xa denotes the possible configurations of subsystems a = A,B, AB). The correla- 
tions between A and B are given by 

I{A: B) = H{A) + H{B) - H{AB). (55) 

We are now ready to formulate and prove the following theorem | [94) . 

Theorem 14. Let L be a lattice with d-dimensional classical spins at each site. Let 
p be a Gibbs probability distribution coming from finite-range interactions on L. 
Then, dividing L into regions A and B, one has 

I{A: B) <\dA\\ogd. (56) 

Proof. First, notice that the Gibbs distributions coming from finite-range interac- 
tions have the property that if a region C separates A from B in the sense that no 
interaction is between A and B then p{xa\xc, xb) = p{xa\xc)7 which we rewrite 
as 
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p{xA,xc)p{xB,Xc) 

P[Xa,xb,xc)^ 7 — ^ • (57) 

P[xcj 

Now, let A and B be two regions of L, and let dA C A and dB C Bhe boundaries 
of A and B, respectively, collecting all sites interacting with their exteriors. Finally, 
let A = A \ dA and B = B \ dB. Since OA separates A from dB (there is no 
interaction between A and dB), we can use Eq. ( 57 1 to obtain 

H{AB) = H{AdAB) = - ^ p{xA,XQA,XB)\og2p{xA,X0A,XB) 

= - X! P(2^A7a;aA)log2P(x^,a;aA) 
^ X! P(2^9^'2;B)log2?3(a;a^,a;B) 

+ X! ^(2^9-4) l0g2 ^(2:9^) 

xaA 

^H{A) + H(dAB)-H(dA). (58) 

Since 9i3 separates OA from i?, the same reasoning may be applied to the sec- 
ond term of the right-hand side of the above, obtaining H{dAB) = H{dAdB) + 
H{B) - H{dB). This, together with Eq. (|58]), gives 

H{AB) = H{A) + H{B) + H{dAdB) - H{dA) - H{dB), (59) 

which in turn after application to Eq. ( [55] l allows us to write 

I{A: B) ^ I{dA: dB). (60) 

It means that whenever the probability distribution p has the above Markov property, 
correlations between A and B are the same as between their boundaries. 

Now, we know that the mutual information can be expressed through the con- 
ditional Shannon entrop>[^as I{X : Y) = H{X) - HiX\Y). Since H{X\Y) is 
always nonnegative, we have the following inequality 

I{dA : dB) < H{dA) log d. (61) 



To get Eq. (56i it suffices to notice that H{A) is upper bounded by the Shannon 
entropy of independently and identically distributed probability p{xa) — l/d)^^, 
which means that H{A) < \A\ log d. ■ 

Let us now show that a similar conclusion can be drawn in the case of quantum 
thermal states |94|, where the Markov property does not hold in general. 

Theorem 15. Let L be a lattice consisting of d-dimensional quantum systems di- 
vided into parts A and B (L — AU B). Thermal states (T > 0) of local Hamilto- 



The conditional Shannon entropy is defined as H(A\B) = H{A, B) — H{B). 
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nians H obey the following area law 

I{A:B)<l3iT[Ha{QA®QB~QAB)]- (62) 

Proof. The thermal state gp = e~^^ /tv{e^^^) minimizes the free energy F{g) = 
tT{Hg) - {l/l3)S{g), and therefore F{gp) < F{gj ® g^) with g^ and g^ subsys- 
tems of gi3. This allows us to estimate the entropy of the thermal state as 

Sigfi)^p[tTiHgp)~F{gp)] 

> P [ir{Hgp) - F{gj ® g^^)] 

= /? [ir{Hgp) - triHgf (g> gj^)] + S{g^ ® pf) 

= /? [tr{Hgp) - tr{Hgj (g> g^)] + Sigj) + Sig^), (63) 

where the last equality follows from additivity of the von Neumann entropy S{p ^ 
a) — S{p) + S{a). Putting Eq. ([63]l into the formula for mutual information we get 



I{A ■.B)<(3 [triHg^ ® g^) - triHg,,)] . (64) 

Let us now write the Hamiltonian as H = Ha + Hq + Hb, where Ha and Hg 
denote all the interaction terms within the regions A and B, respectively, while Hq 
stands for interaction terms connecting these two regions. Then one immediately 
notices that tr[iJ^( b){q'b ® — Qp)] =0 and only the Hg part of the Hamiltonian 



H contributes to the right-hand side of Eq. (64 1. This finishes the proof 



Let us notice that the right-hand side of Eq. (62 1 depends only on the boundary, and 



therefore it gives a scaling of mutual information similar to the classical case ( [6T] ) 



Moreover, for the nearest-neighbor interaction, Eq. (62i simplifies to I{A : B) < 
2/3 11 ft, 1 1 \dA\ with \\h\\ denoting the largest eigenvalue of all terms of H crossing the 
boundary. 



4 The tensor network product world 

Quantum many-body systems are, in general, difficult to describe: specifying an 
arbitrary state of a system with A^-two level subsystems requires 2^ complex num- 
bers. For a classical computer, this presents not only storage problems, but also 
computational ones, since simple operations like calculating the expectation value 
of an observable would require an exponential number of operations. However, we 
know that completely separable states can be described with about N parameters — 
indeed, they correspond to classical states. Therefore, what makes a quantum state 
difficult to describe are quantum correlations, or entanglement. We saw already that 
even if in general the entropy of a subsystem of an arbitrary state is proportional 
to the volume, there are some special states which obey an entropic area law. In- 
tuitively, and given the close relation between entropy and information, we could 
expect that states that follow an area law can be described (at least approximately) 
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with much less information than a general state. We also know that such low entan- 
glement states are few, albeit interesting — we only need an efficient and practical 
way to describe and parameterize them. 



4.1 The tensor network representation of quantum states 

Consider a general state of a system with N d-level particles, 

d 

IV') = X! Ciii2...»ivl*li*2, • ■ • ,«Af)- (65) 
11,12, ■■■,««=! 

When the state has no entanglement, then Ci.i, i., = c,-^''c,-^^ . . . c,-^"* where all c's 
are scalars. The locality of the information (the set of coefficients c for each site is 
independent of the others) is key to the efficiency with which separable states can 
be represented. How can we keep this locality while adding complexity to the state, 
possibly in the form of correlations but only to nearest-neighbors? As we shall see, 
we can do this by using a tensor at each site of our lattice, with one index of the 
tensor for every physical neighbor of the site, and another index for the physical 
states of the particle. For example, in a one-dimensional chain we would assign a 
matrix for each state of each particle, and the full quantum state would write as 



tr 



|ll,i2,...lAr), (66) 



where stands for a matrix with dimensions Dk x Dk+i- A useful way of un- 
derstanding the motivations for this representation is to think of a valence bond pic- 



ture 195) . Imagine that we replace every particle at the lattice by a pair (or more in 
higher dimensions) of particles of dimensions D that are in a maximally entangled 
state with their corresponding partners in a neighboring site (see Figure]?]). Then, 
by applying a map from this virtual particles into the real ones, 

d D 

^ = E E (67) 

i=l a, 13=1 



we obtain a state that is expressed as Eq. (66 1. One can show that any state |?/') e 
can be written in this way with D = max„i I?„i < d^/^. Furthermore, a 
matrix product state can always be found such that |96j 



Af'^^Af^ = Id., for 1 < fc < N, 

A|['=U['=-iUfJ = ylW, for 1 < fc < iV, and 
For open boundary conditions yl^^l = yl^^l 1, and yll'^'l is a Z3fc+i x Dfe+i 
positive diagonal matrix, full rank, with tvA^^^ = 1. 
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Fig. 4 Schematic representation of tensor networks. In panel (a) we stiow tlie meaning of the 
elements in the representation, namely the solid line joining two virtual particles in different sites 
means the maximally entangled state between them, and the grey circle represents the map from 
virtual particles in the same site to the physical index. In panel (b) we see a one-dimensional tensor 
network or MPS, while in (c) we show how the scheme can be extended intuitively to higher 
dimensions — in the two-dimensional example shown here, a PEPS that contains four virtual 
particles per physical site. 



In fact, yll'^'l is a matrix whose diagonal components Af^, n = 1, . . . , Dk, are the non- 
zero eigenvalues of the reduced density matrix obtained by tracing out the particles 
from fc + 1 to A^, i.e., the Schmidt coefficients of a bipartition of the system at site 
k. A MPS with these properties is said to be in its canonical form [91] . 



Therefore, Eq. (66 1 is a representation of all possible states — still cumbersome. 
It becomes an efficient representation when the virtual bond dimension D is small, 
in which case it is typically said that the state has a matrix product state (MPS) 
representation. In higher dimensions we talk about projected entangled pair states 
(PEPS) |98|. When entanglement is small (but finite), most of the Schmidt coef- 
ficients are either zero or decay rapidly to zero |82J. Then, if contains little 
entanglement, we can obtain a very good approximation to it by truncating the ma- 
trices j4 to a rank D much smaller than the maximum allowed by the above theorem, 
In fact, we can demonstrate the following 



Lemma 1. |97j/ There exists a MPS IV'd) with bond dimension D such that || ji/)) — 
WdW < 2EL7 ^»iD), where e^{D) = eS+T""' ' 

Proof. Let us assume that the MPS is in its canonical form with D = 2^^^. Defining 
a projector into the virtual bond dimension P = EfcLi ^ TPCM map 

KiX) = J2^ vlf'^^XAf"!, we can write the overlap 



{lb\lbD) = Tr [$2(. . . $N-2{$N-iiA^''-'^P)P)P) . . .)P 



(68) 
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By defining = $fe(y['=+ilp), with yl^-i) = yll^-^lp, and using that 

Tr|$(X)| < Tt\X\, we can see that 

tr|ylW -yW| tr|$fe(yl['=+il -yl'^+ilp)! 

< tr|yl[''^+il +tr|yl['^-+il(l-P)|, (69) 

where the last term is equal to 22a=D+i ■ Finally^ applying this last inequality 
recursively from — 1 to 2, and using that (V'd IV'-d) ^ !> can obtain the desired 
bound on {ipltpo). □ 

Lemma ([T]i is most powerful in the context of numerical simulations of quantum 
states: it gives a controllable handle on the precision of the approximation by MPS. 
In practical terms, for the representation to be efficient the Schmidt coefficients A 
need to decay faster than polynomially. However, we can be more precise and give 
bounds on the error of the approximation in terms of entropies l99] : 

Lemma 2. Let Sa{p) = log(tr/3")/(l — a) be the Renyi entropy of a reduced den- 
sity matrix p, with < a < 1. Denote e{D) = X^i^D+i ^i' ^'th Aj being the 
eigenvalues of p in nonincreasing order Then, 

\og{e{D)) < — (s^ip) - log . (70) 

The question now is when can we find systems with relevant states that can be 
written efficiently as a MPS; i.e. how broad is the simulability of quantum states by 
MPS. For example, one case of interest where we could expect the method to fail is 
near quantum critical points where correlations (and entanglement) are singular and 
might diverge. However, at least in ID systems, we can state the following: 

Lemma 3. [97] In one dimension there exists a scalable, efficient MPS representa- 
tion of ground states even at criticality. 

Proof. In one dimension, the worst case growth of entropy of a subsystem of size 
L, exactly at criticality, is given by 

5„(pi)^^(^l + i^logL. (71) 

Let us take the length L to be half the chain, N = 2L. By means of the previous 
discussion, we can find a MPS \iPd) such that its distance with the ground state 
\tpGs) is bounded as || \iPgs)~'\''I'd) IP < co/^' with eq constant. Now, let be the 
minimal virtual bond dimension needed for this precision, i.e. || \ipGs) ~ iV'z?) IP 5; 
2 X 2L e„iax{D). We demand that 



f-max{D) < 

< exp 



log L — log 



12 a 1-a 



(72) 
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Set 


const 


logL 


L''{k < 1) 


L 


Sc. < 1 


OK 


OK 


9 


7 


5 = Si 


7 


? 


7 


NO 


Sc. > 1 


? 


? 


NO 


NO 



Table 1 Relation between scaling of block Renyi entropies and approximability by MPS |99| ). In 
the "undetermined" region denoted with question marks, nothing can be said about approximability 
just from looking at the scaling. 



from which we can extract 

< const ^^^^ LTi^^ cx poly(L). (73) 



Establishing that there exists an efficient representation of the ground state is 
not enough: we must also know if it is possible to find it efficiently too. In one 
dimensional gapped systems, the gap A typically scales polynomially, which means 
that DMRG and MPS methods should converge reasonably fast. One can, however, 
formalize the regime of efficiency of MPS as a function of how the different Renyi 
entropies scale with subsystem size | [99) . In table ([T]i we summarize the currently 
known regimes where the MPS approach is an appropriate one or not. 



4.2 Examples 

Here we present a few models and states that are fine examples of the power of MPS 
representations |97|. 

Example 6. A well known model with a finite excitation gap and exponentially de- 
caying spin correlation functions was introduced by Affleck, Kennedy, Lieb, and 
Tasaki | |100||10l| — the so called AKLT model. The model Hamiltonian is 

if - ^ • + ^ {S, ■ S,+if . (74) 

i 

For 5 = 1, the local Hilbert space of each spin has three states, thus d = 3. The 
ground state of this Hamiltonian can be written compactly using a translationally 
invariant MPS with bond dimension D = 2, specifically 

A_i = <T^, Aq = V2a+, Ai = -V2a~. (75) 



Example 7. A paradigmatic example of a frustrated one dimensional spin chain is 
the Majumdar-Ghosh 1 102| model with nearest and next nearest-neighbor interac- 
tions: 
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-ff = ^ 2(7, • tTj+i + (Ti ■ (Ti+2, (76) 

i 

The model is equivalent to the or Ji — J2 Heisenberg model with J1/J2 = 2. 
The ground state of this model is composed of singlets between nearest-neighbor 
spins. However, since the state must be translationally invariant, we must include 
a superposition of singlets between even-odd spins, and "shifted" singlets between 
odd-even spins. The state can be written compactly in MPS form using D = 3, 



^0 = -1 , = 1 (77) 





Example 8. A relevant state for quantum information theory is the Greenberger- 
Home-ZeiUnger (GHZ) state, which for N spin 1/2 particles can be written as 

10)®^ + \1)<»N 

\GHZ) = ^ . (78) 

v2 

GHZ states are considered important because for many entanglement measures they 
are maximally entangled, however by measuring or tracing out any qubit a classical 
state is obtained (although with correlations). GHZ states can be written using D = 
2 MPS, specifically Aq^i = 1 ± a^- Also the "antiferromagnetic" GHZ state is 
simple, Aq^i = u^. 

Example 9. Cluster states are relevant for one-way quantum computing. They are 
the ground state of 

// = ^af_i<7fc7|+i, (79) 
and can be represented using &D = 2 MPS, 

^=i;;).-^-(i-c') 

Example 10. (Classical superposition MPS) Imagine we have a classical Hamilto- 
nian 

iJ = ^/i((Ji,a,), (81) 

where ai = 1, . . . , d, and h{ai,aj) are local interactions. The partition fimction of 
such a model at a given inverse temperature ^ is 

Z = J2exp[-pH{a)], (82) 

where the sum is over all possible configurations of the vector a. Let us now define 
a quantum state whose amphtude for a given state of the computational basis 
corresponds to the term in the partition function for that state, i.e. 
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|(7l . . . CTat) 



13 



h{ai,aj) 



|(7i . . .CTat). 



(83) 



We shall now define a map P — in the same manner as in valence bond states — that 
goes from to such that 



P\s,k) = \s){(ps\k), 



where we have defined 



^{iPs\a){ips\a) = exp 



a=l 



(84) 



(85) 



To visuahze what happens when we insert these back into the classical superposition 
state li'/j), let us concentrate for a moment on a one-dimensional system: 



1 

71, 



E ®^p 



/3 



h{ai,a2) 



. . . exp 



-/l((Tjv-i,crjv) 



\ai . . .Gn), 



era 1*^2/ 

(Ti,.")<7Af \ai — 1 a2 = l 

d \ 

X ... X E (i^<^^|ajv)(v'CTi|Q!iv) I ki . . .crjv), 

a;iv = l / 

1 

= ^ E E [(V<^il"Jv)(</'<^il"l)] [(V'<72|Q!l)(</5a2|a2)] 



(86) 



X ... X [((/?cr„|aAr-l)(</'<Tjv|aA')] ki ■ ■ ■<TAr),(87) 



and we can replace ^g*^„ ^ = (i^Si|Q)(</'sJ/3), thus expressing the classical thermal 
superposition state as a MPS. 

These states have some important properties: 

(i) They obey strict area laws, 

(ii) They allow to calculate classical and quantum correlations, and 

(iii) They are ground states of local Hamiltonians. 

Property (i) should be obvious by now, since we have explicitly shown the MPS 
form of the state. We can show easily property (ii) for Ising models. A classical 
correlation function / must be evaluated with the partition function, (/((t)) = 
Yl,a /(o')e~'^^^'^) jZ, but this is just the expectation value of an operator made of 
changing the argument of / into operators, and evaluated with |V'/3)- Since it is 
the expectation value of a MPS, it is efficient to compute. Finally, we will demon- 
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strate property (iii) at length in the next section, because it will lead us into the final 
topic of this lectures: Quantum kinetic models. 



4.3 Classical kinetic models 



Our goal in this section is to show the local Hamiltonians whose ground state is 
the classical superposition state defined in the previous section. As we shall see, 
these Hamiltonians will arise from the master equation of a classical system that is 
interesting in its own right, so we will first spend some time on it. 

Let us consider a system made out of N classical spins interacting through a 
Hamiltonian H. If cTj denotes the state of ith spin, we will label the configurations of 
the system by cr = (tri, . . . , ctat), and the probability of finding at time t the system 
in state a (given that it was in state ctq at time to) by P{<j, t) = P{a, t\ao, to). In 
what follows we focus on this probability distribution, whose dynamics is described 
by a master equation: 

cr' cr^ 

where W((T, a') is the transition probability from state a' to state a. This equation 
defines the class of kinetic models, and it clearly describes a Markov process — the 
instantaneous change of P{a, t) does not depend on its history. 

We will only consider systems that obey a detailed balance condition, i.e. 

VF(a,a')e-'3^("') = W{a\a)e-l^"^^\ (89) 

With this condition, the stationary state of the master equation (the one that fulfills 
Pst{<J^ t) — 0) is simply Pst{<y) — e^^^^"'' jZ, with Z being the partition function. 
This state in particular will map into the classical superposition state defined above. 



but we still have not found its parent Hamiltonian. For this, we will rewrite Eq. ( 88 i 
in the form of a matrix Schrodinger equation (albeit with imaginary time) from 
which we can identify a Hamiltonian. 

Let us apply the transformation ■(/'(c, t) — e^^(°'^/^P((T, <), which leads to 

i) - E e^^(")/'VK((T, (T')e-'^^("')/2^(a', t) - W{g' , a)^{<j, t) 

a' 

= -J2H0ia,a')i;ia\t) (90) 

a' 

with 
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Notice that the detailed balance condition guarantees that the matrix Hj3 is Hermi- 
tian, so we can interpret it as a Hamiltonian. Furthermore, because of the conserva- 
tion of probability. Hp can only have non-negative eigenvalues, which means that 
the state associated to the stationary state Pgt with eigenvalue zero must be a 
ground state: the classical superposition MPS that we were looking for 

Remarkably, we have not said anything yet about H. A famous example of such 
kinetic model is a single spin-flip model considered by Glauber | ,103J , for which H 
is the Ising Hamiltoniarj^ 

Hia) = ifising(a) ^-JJ2 ^1 (-^ > 0)- (92) 
Denoting by Pi the flip operator of the i-th spin, i.e., f icr, ~ —(Ji, the general master 



equation ( 88 1 reduces in this case to 

P((T, t) = J2 [Wia, P^a)PiP,a, t) ~ W{P,a, a)P{a, t)] (93) 

i 

with W{a,Pia) now called spin rates. It was shown in flOSl that the most gen- 
eral form of spin rates with symmetric interaction with both nearest-neighbors, and 
satisfying the detailed balance condition ([89|, is given by 



w{P,cj,a) = r(l + - (l/2)7a,(^,_i + (94) 

with r > 0, —1 < 5 < 1, and < 7 < 1- The S — case was thoroughly 
investigated by Glauber |103| , who showed that all the relevant quantities can be 
derived analytically — including the dynamical exponent that turned out to be z = 



2. The more general case of nonzero S was treated in a series of papers 1 104- 106 1, 
that showed for instance that the choice 6 = 7/(2 — 7) leads to an interesting 
dynamical exponent z 7^ 2. 

If we rewrite the single spin-flip master equation in the form of the Schrodinger 
equation, we obtain an associated quantum Hamiltonian 



77^(5,7) ^-rJ2 - S(<5,7)<-i<+i) < 



-(1 + daU^t+^ (1 - (l/2)7^f i^Li + <Tf+i))] , 

(95) 



where 



A(J,7)-— ^^^;^=r-^, B{S,j)^l^A{S,j) (96) 
2(1- VI -7^) 

and CT^ and are the standard Pauli matrices. For 6 = this Hamiltonian was 
diagonalized in Ref. |107) , and independently in Ref. ]108| . 



' This is the reason why the Glauber model is also known as the kinetic Ising model (KIM). 
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The Hamiltonian 7), and also the other ones that can be derived in this 

way, are typically gapped except at a critical temperature where the gap vanishes 
with the critical exponent z that characterizes the model. In one dimension /3c ~ 00, 
but for larger dimensions this model has a finite critical temperature. 

We have seen thus far how the the master equation of a classical spin model (that 
obeys the detailed balance condition) can be associated to a quantum Hamiltonian 
with some interesting critical properties — for example, its ground state obeys a 
strict area law and can be written efficiently as a MPS. Nevertheless, the underlying 
model is still classical. In the next section, we will see one way in which we can 
generalize the initial model to be quantum, while retaining the same structure that 
leads to associated Hamiltonians that obey area laws. 



5 Quantum kinetic Ising models 

Here we discuss ways to generalize the kinetic equation ( [88] l to a quantum master 
equation, but in such a way that its diagonal part reproduces the corresponding ki- 
netic model. A similar approach was taken in Ref. |jT09|, where a quantum master 
equation that reproduced a kinetic Ising model was proposed (see also Ref. |110| ). 
However, no attempts aiming at fully solving such QMEs are known so far Our 
purpose is to give quantum generalizations of the classical kinetic models that can 
be solved analytically. 

Recently, we presented such a generalization |1 1 Ij fo r the single spin-flip model, 
Eq. ( 93 I, with the spin rates of Eq. ( 94 1. In Ref. jl ll") we were able to decouple 



the master equation for the density matrix of a quantum system into 2 
equations with the same structure as the ones studied above. Here, we will only show 
the associated Hamiltonians (and their spectra) obtained in these models. However, 
we will demonstrate how to approach the problem but in a different model that 
allows transitions that flip two consecutive spins. 



5.1 A two spin flip model 



First, let us particularize the classical kinetic equation (88 1 to the case where the flip 



operator acts on pairs of consecutive spins of the chain, i.e. 

(97) 

where Fi^i+i denote spin flips at positions i and i + 1, while the spin rates are 
given by Wi{Fi^i+ia,cr) = r[l - {l/2)-f{ai^iai + ai+iai+2)] with < T < 00 



and 7 = tanh2/3J. This model was investigated in Ref. ]1 12| , where the associ- 
ated Hamiltonian was found and diagonalized using the Jordan-Wigner transforma- 
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tion |1 13| followed by Fourier and Bogoliubov-Valatin |114||115| transformations. 
In particular, Hilhorst et al. were able to show that, despite the complexity of the 
transformations, one can easily compute expectation values such as magnetization, 
energy density, or correlations, and that they have a relatively simple exponential 
behavior 01121. 



Here we will define through a master equation a quantum model that resembles 
the kinetic model above. For this, we will replace classical probabilities with the 
quantum density matrix, and classical operators with quantum ones (e.g. cr^ is the 
qubit flip operator). Consider the following master equation 



dtQ{t) = J2 



(98) 

where {•, •} denotes the anticommutator and Wi{a^) are quantum mechanical gen- 



eralizations of the spin rates ( 94 1, now written in terms of the operators. 



w.{an ^r[l- ^(^ti^f + <^f+i<2)] • (99) 

Although it looks complicated, the quantum kinetic model above can still be 
solved with techniques similar to the classical case 1111]: the key ingredient is to 
find a large number of constants of motion that allow to split the master equation into 
a set of ordinary Schrodinger equations. To see this, we must represent the density 
matrix g{t) as a vector in an expanded Hilbert space. This follows from a simple 
isomorphism between linear operators from AId{C) and vectors from . In other 
words, writing our density matrix in the computational basis in (C^)®^ as g{t) — 
Ea.sbWl-T.skXc^l' we can treat it as a vector \g{t)) = E^,;? bWlrr.slcr) |(?> from 
(C^)®^ (g) (C^)®^. Even if formally we are enlarging the number of spins from 
N to 27V, the advantage is that now we deal with "pure states" instead of density 
matrices which allows us to find many conserved quantities. This, in turn, shows 
that the effective Hilbert space used is much smaller than the initial one. To be 
consistent, operators that appear to the right of g{t) must be replaced with "tilded" 
operators that act on the right subsystem of the expanded space, while operators on 
the left of the density matrix ("untilded") act on the left subsystem (for instance 
^fc^f I'S)!^ = This notation allows us to rewrite the master equation 



( 98 1 as the following matrix equation 

(100) 

As was the case for the initial classical master equation, the matrix appearing on the 



right-hand side of Eq. (100 1 is not Hermitian. In order to bring it to Hermitian form 



we can use the detailed balance condition, which suggests the transformation 

|£)(t)) =cxp[-(/3/4)[H((T)+?{(5)]]|VX0), (101) 
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with y, denoting the quantum generahzation of the Ising Hamiltonian H — —JJ2i '^i'^i+i- 
With this transformation, and denoting 



(102) 



Eq. ( 100 1 can be written as 



\m) = E 



\m) 



(103) 

which we can see as a Schrodinger equation \ilj{t)) = —H\ip{t)) with Hermitian H. 

We have reached the point where all these changes of notation payoff: indeed, 
the form of H makes it clear that it commutes with af afj^iafafj^i (i = 1, . . . , N). 
Therefore, we can introduce new variables = '''f ^'f+i^^f 2r,f_^i (* — 1, ■ • • , ^) 
which are constants of motion and reduce the number of degrees of freedom. 
In particular, tilded variables can be expressed by a and the new variables r as 
^i^i+i — foJ" ™y 111 Other words, we have replaced a and cr by r and a, 

of which r is conserved. To each configuration of r's we associate a natural num- 
ber from to 2^ — 1, which corresponds to a particular correlation between the 
a and a variables. For example, t = corresponds to all r-spins up (r^ = 1 for 
i = 1,...,N), while r = 2^ - 1 means that n = -1 for i = 1, . . . , iV. With 
this notation, each value of r is associated to a Hamiltonian Hr that acts only in the 
space of N spins and is of the form 



E 



1, 



(104) 



where rcr^ denotes r^CTf {i — 1, . . . , N). Because these Hamiltonians are indepen- 
dent from each other, we have converted the problem of solving the general master 
equation ( 98 i to the problem of diagonalizing 2^ Hamiltonians, each of dimension 
2^ X 2^. Now, we have that 



\Mt)), 



2"-l 



(105) 



After simple algebra one sees that the explicit form of Hr is 



, (106) 



where 
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COS^Lp, 



with 



cos (p ■ 



n-in+i = 1 

^/cOs2^p, Tj-in+i = -1 



cosh/3J 
(cosh2/3J + smh2/3J)i/2' 



sin ^p, n^iTi+i = 1 



0, 



Ti-in+i = — 1, 

(107) 



sm = 



sinh/3 J 



(cosh^/?J + sinh^/?J)i/2' 



(108) 

and f{x) = (1/2)(1 + a;). Here the angle ranges from zero (which corresponds 
to infinite temperature) to 7r/4 (which corresponds to T = 0) and in this notation 
7 = sin2(^. Let us notice that for r = Eqs. (106i and (107i, as it should be, 
reproduce the Hamiltonian derived in |1 12|. This, however, contrary to the single 
spin-flip case, is not the case for r = 2^ — 1, where one of the terms in the square 
brackets vanishes and the Hamiltonian reduces to 



i+2l 



'i+l 



-1 



(109) 



Let us discuss now some of the properties of Hr- Below we show that for all r 
they are always positive operators. We also find all the cases with respect to if and 
T for which the Hamiltonians can have zero-energy ground states. 

Lemma 4. The Hamiltonians Hr are positive for any t = 0, . . . , 2^ — 1. 



106 1. The 



Proof. Let us denote by H)- ' the ith term appearing in the sum in Eq. 
idea is to show that all H)- are positive, and the positivity of Ht follows immedi- 
ately. Of course, the form of i?!''' changes depending on r-spins at positions i — 1 
and i + l. Therefore, we distinguish several cases with respect to different possible 
configurations of these spins. 



For ri_i — Tj+i, one easily infers from Eqs. ( 106 1 and ( 107 1 that 



- (cosV - sinV'Tf^irTfcrf+iCrf+a) (HO) 

In the case when both spins Tj_i and r^+i are down, the function / is zero and 
both terms in square brackets vanish and the above operator becomes 1 — (cos^ </? — 
sin^ <^i^i<^i(^i+i<^i+2)'^i '^i+i- 1'^ Is clear then that its minimal eigenvalue is zero. 
In the case when Ti_i = t^+i = 1 these terms do not vanish, however, still this is 
effectively a 16 x 16 matrix which can be shown to be positive computationally: 
using the software Mathematica we can easily see that the minimal eigenvalue is 
zero. 

For Ti_i = — Ti+i, one of the values f{Ti^i) or /(ri+i) is zero. Assuming that 
f{Ti-i) — (the case of /(ri+i) = leads to the same eigenvalues), one has 

= 1 - ^ sin 2^ - v/^^ <af+i . (Ill) 
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(i) 

When constrained to three consecutive spins (i — 1, i, and i + this Hr ' is just a 8 
by 8 matrix (on the remaining spins it acts as the identity matrix) and its eigenvalues 
can be obtained using Mathematica. One then checks that its minimal eigenvalue is 
1 — (l/2)-\/4cos2(y9 + (sin2(y9)2 with ip e [0, 7r/4]. Simple analysis shows that this 
is a nonnegative function of ip and gives zero only when ip = 0. 

In conclusion, Hr'' > for all rs and (p e [0, 7r/4] and therefore our Hamiltoni- 
ans Hr are positive. ■ 

Based on the above analysis, let us now distinguish all the cases with respect to 
r and (p when Hr{f) can have zero-energy eigenstates. It clearly follows from the 
proof of lemma|4]that if r 7^ or r 7^ 2^ — 1 there exists i such that rj_i ^ r^+i 
and then the corresponding Hr{<p) can have zero eigenvalues only when ip = Q. Let 



us now discuss this case. It follows from Eqs. ( 106 1 and ( 107 1 that for (,5 = 0, which 



corresponds to infinite temperature, the dependence on r vanishes and one obtains 
7J.(0)^iT = ^(l (112) 

i 

which has a doubly degenerate ferromagnetic ground state. 
For r = one gets the Hamiltonian obtained in |1 12j , that is 

i 

- (1 - (l/2)sin2^(af„iaf . (113) 

The ground state of this Hamiltonian is doubly degenerate for all values of p, ex- 
cept for ip — tt/A (zero temperature) where also the first excited state becomes 
degenerate with the ground state [lTT\. For many values of r this statement holds, 
except that the ground state has a positive energy — implying that the off-diagonal 
elements of the QME decay in time. In other cases we find that the ground state is 
unique for all values of p, even tt/4. Typical spectra for some values of t in finite 
systems are shown in Fig. [5] 



5.2 The single flip model 

For comparison only, we reproduce here the associated Hamiltonians that are ob- 
tained when single flip processes are allowed in the quantum master equation [1 1 1[ . 
Again, a set of conserved quantities allows us to break the QME into 2^ Schrodinger 
equations labeled by a parameter r 

\Mt)) = -HrlMt)) (t = G,...,2^-1), (114) 



where the Hamiltonians Hr are given by 
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Fig. 5 Low energy states of the Hamiltonians dl06k associated to the two flip quantum master 
equation for a system with N = Id spins as a function of tp. The panels are (a) r = 2* — 1 (half 
T-spins up and half down), (b) r = 2* (only one r-spin up, the others down), and (c) r = 2^ + 2* 
(two neighboring r-spins up, the others down). Only in case (c) the ground state is fully degenerate 
for all values of Lp, in the other two the first excited state energy is very close but not equal to the 
ground state. In case (b) the ground state is not degenerate at (/p = 7r/4, while in the other two cases 
it is. 



i 

-(5/(T,_iTi+i)CTf_iCTf+i] , (115) 



where 
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(116) 



(117) 



Here, each r means a configuration of the conserved quantities that is different than 
the one shown above for the two flip model - however, we still use its binary repre- 
sentation so that T is a shorthand notation for N variables {ti, . . . ,tn), each taking 



values ±1. The equation in ( 114 1 for r = corresponds to the diagonal elements of 
g{t), while for the remaining t 7^ 0, they describe the off-diagonal elements of the 
density matrix. 

Let us shortly comment on th e abo v e mo del. F irst, i t is easy to notice that for 
T = or T = 2^ — 1, from Eqs. ( 1 15 1, ( 1 16 1, and ( 1 17 1 one recovers the Hamilto- 
nian ( 95 1. Since, as shown in Ref. 1 107 1 the Hamiltonian ( 95 1 has a ground state with 



zero energy, it means that there exist off-diagonal elements surviving the evolution. 
On the other hand for r ^ 0, 2^ — 1 one gets (95 1, however, with some impurities. 
After substitution of bond variables (see e.g. 1 108 1) one can map to disordered 
Heisenberg chains meaning that for some particular values of the involved parame- 
ters the model can be solved analytically. On the other hand, one may always treat 
this model numerically through matrix product states. 

We show in Fig. |6]the spectra for some of these Hamiltonians, which is to be 
contrasted with the spectra from the two spin flip models. Fig. [5] In the single flip 
model the ground state is always unique except at zero temperature, where for all of 
the associated Hamiltonians one observes criticality. 



6 Discussion and Outlook 

In these lectures we have seen how quantum information theory can bring about a 
fresh perspective into many-body physics. However, the field is much bigger than 
what we have reviewed. Let us just mention here a few relevant topics that we have 
not covered, but that have received plenty of attention from the community, and that 
certainly have contributed to sizable advances in our understanding of many body 
physics. 

One interesting application of entanglement is to critical phenomena. We briefly 
saw how block entanglement entropy scales differently at a gapless critical point. 
However, many entanglement measures display some kind of special behavior 
around quantum criticality — which was first observed Q in the concurrence of 
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Fig. 6 Low energy states of the single flip Hamiltonians studied in Ref. as a function of the 
temperature parameter 7 = tanh2/3J for the same parameters as in figure [5] The three panels 
correspond to the same r-spin configurations, even though the variables r are defined differently. 
Notice that in this case the spectra becomes degenerate always at 7 = 1, and that the ground state 
is always unique. 



nearest-neighbor spins of an Ising chain (see Ref. | |116| for a recent review of ac- 
tivity in this field). Quantum criticaUty, in fact, is a very active subject in the con- 
densed matter community, and has been studied using other quantum information 
approaches hke the ground state fidelity 1 117| and the Loschmidt echo 1 1 18 1, whose 
usefulness in practice has been demonstrated experimentally 1 ,1 19 120] . 

Another problem that is gaining interest is that of topological order, which we 
mentioned briefly as one of the motivations for studying area laws. One interesting 
recent development is the study of entanglement spectra p2ip22) , defined through 
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the Schmidt decomposition in such a way that each Schmidt coefficient of a 
bipaitition is interpreted as a dimensionless energy = — log A^. This approach 
allows to generalize the von Neumann block entropy by introducing a virtual tem- 
perature, and study the structure of entanglement with more detail. In particular, it 
appears that gapped systems with topological order always have a gapless entangle- 
ment spectrum |121|. As a characterization of entanglement, the whole spectrum 
promises to be better than just entropy — simply because a set of numbers contains 
more information than a single one. 

Although we concentrated mostly on the theoretical aspect of matrix and ten- 
sor product states, the field is also strongly geared to the practical application of 
simulation of many-body systems in classical computers. On the theory side, the 
tensor product approach has given successful advances in the theory of computa- 
tional complexity applied to quantum mechanics |99|, and in the recent theory of 
entanglement renormalization 1 123). 



On the computational side, MPS algorithms have expanded the effective DMRG 
methods, and tremendous progress is being done in the simulation of strongly corre- 
lated particle systems. Bosonic particles can be represented straightforwardly |124) 
by mapping the d internal levels of the spins to the occupation number at each lat- 
tice site - thus truncating the Hilbert space to the subspace with at most d — 1 
particles in each site. Fermionic models, however, require some extra care when 
contracting the indices in the network so that fermionic commutation relations are 
respected | 125f |130|. In any case, tensor networks are still computationally efficient 
with strongly correlated electron systems, which puts these algorithms at an ad- 
vantage over quantum Montecarlo type techniques - who suffer from the so called 
"sign problem" in this type of systems 1 131 1. Therefore, tensor network techniques 
are important to support the large experimental efforts towards implementing quan- 
tum simulations of fermionic models (mainly with trapped ions y_32j and ultracold 
atomic systems |133| ). 

The quantum kinetic Ising models discussed in the last section hold plenty of 
potential for the near future in at least two fronts. First, they represent a whole new 
class of many-body systems amenable to analytical solution, and can therefore bring 
new insight into our understanding of complex quantum many-body dynamics, as 
well as some classical reaction-diffusion problems | ,134J . Second, they have a close 
relationship and could be useful to the recent ideas on "environment design" ]135j - 
|137) : crafting and/or manipulating the environment of a system so that it is driven 
to an interesting quantum many-body state, usually with a dynamics given by a 
quantum master equation. Because quantum kinetic models can be well understood 
and controlled, they might provide the foundation on top of which more elaborated 
systems are designed. 



Acknowledgements We are grateful to LI. Masanes for helpful discussion. We acknowledge the 
support of Spanish MEC/MINCIN projects TOQATA (FIS2008-00784) and QOIT (Consolider 
Ingenio 2010), ESF/MEC project FERMIX (FIS2007-29996-E), EU Integrated Project SCALA, 
EU STREP project NAMEQUAM, ERC Advanced Grant QUAGATUA, Caixa Manresa, AQUTE, 
and Alexander von Humboldt Foundation Senior Research Prize. 



46 



R. Augusiak, F. M. Cucchietti, and M. Lewenstein 



References 

1. D. Jaksch, H.-J. Briegel, J. I. Cirac, C. W. Gardiner, and P. ZoUer, Entanglement of Atoms 
via Cold Controlled Collisions, Phys. Rev. Lett. 82, 1975 (1999). 

2. T. J. Osborne and M. A. Nielsen, Entanglement, Quantum Phase Transitions, and Density 
Matrix Renormalization, Quantum Inf. Proc. 1, 45 (2002). 

3. T. J. Osborne and M. A. Nielsen, Entanglement in a simple quantum phase transition, Phys. 
Rev. A 66, 032110(2002). 

4. A. Osterloh, L. Amico, G. Falci, and R. Fazio, Scaling of entanglement close to a quantum 
phase transition. Nature 416, 608 (2002). 

5. A. K. Ekert, Quantum cryptography based on Bells theorem, Phys. Rev. Lett. 67, 661 (1991). 

6. C. H. Beimett and S. J. Wiesner, Communication via one- and two-particle operators on 
Einstein-Podolsky-Rosen states, Phys. Rev. Lett. 69, 2881 (1992). 

7. C. H. Bennett, G. Brassard, C. Cr6peau, R. Jozsa, A. Peres, and W. K. Wootters, Teleporting 
an unknown quantum state via dual classical and Einstein-Podolsky-Rosen channels, Phys. 
Rev. Lett. 70, 1895 (1993). 

8. R. Horodecki, M. Horodecki, P. Horodecki, and K. Horodecki, Quantum entanglement. Rev. 
Mod. Phys. 81, 865 (2009). 

9. R. F. Werner, Quantum states with Einstein-Podolsky-Rosen correlations admitting a 
hidden-variable model, Phys. Rev. A 40, 4277 (1989). 

10. P. Horodecki, Separability criterion and inseparable mixed states with positive partial trans- 
position, Phys. Lett. A 232, 333 (1997). 

11. M. Zukowski, A. ZeiUnger, M. A. Horne, and A. K. Ekert, Event-ready-detectors Bell exper- 
iment via entanglement swapping, Phys. Rev. Lett. 71, 4287 (1993). 

12. O. Giihne and G. Toth, Entanglement detection, Phys. Rep. 474, 1 (2009). 

13. A. C. Doherty, P. A. Parrilo, and F. M. Spedalieri, Distinguishing Separable and Entangled 
States, Phys. Rev. Lett. 88, 187904 (2002). 

14. F. Hulpke and D. BruB, A two-way algorithm for the entanglement problem, J. Phys. A: Math. 
Gen. 38, 5573 (2005). 

15. L. Gurvits, Classical complexity and quantum entanglement, STOC 69, 448 (2003). 

16. M.-D. Choi, Positive linear maps, Proc. Symp. Pure Math. 38, 583 (1982). 

17. A. Peres, Separability Criterion for Density Matrices, Phys. Rev. Lett 77, 1413 (1996). 

18. M. Horodecki, P. Horodecki, and R. Horodecki, Mixed-State Entanglement and Distillation: 
Is there a Bound Entanglement in Nature?, Phys. Rev. Lett. 80, 5239 (1998). 

19. K. Zyczkowski, P. Horodecki, A. Sanpera, and M. Lewenstein, Volume of the set of separable 
states, Phys. Rev. A 58, 883 (1998). 

20. D. R DiVincenzo, R W. Shor, J. A. Smolin, B. M. Terhal, and A. V. Thapliyal, Evidence for 
bound entangled states with negative partial transpose, Phys. Rev. A 61, 062312 (2000). 

21. W. Diir, J. I. Cirac, M. Lewenstein, and D. BruB, Distillability and partial transposition in 
bipartite systems, Phys. Rev. A 61, 062313 (2000). 

22. M. Horodecki, P. Horodecki, and R. Horodecki, Separability of mixed states: necessary and 
sufficient conditions, Phys. Lett. A 223, 1 (1996). 

23. E. Bishop and D. Bridges, Constructive Analysis (Springer, Berlin, 1985). 

24. B. M. Terhal, Bell inequalities and the separability criterion, Phys. Lett. A 271, 319 (2000). 

25. V. Gorini, A. Kossakowski, and E. C. G. Sudarshan, Completely positive dynamical semi- 
groups ofN-level systems, J. Math. Phys. 17, 821 (1976). 

26. M.-D. Choi, Completely positive linear maps on complex matrices. Linear Alg. Appl. 10, 
285 (1975). 

27. K. Kraus, States, Effects and Operations: Fundamental Notions of Quantum Theory 
(Springer Verlag, 1983). 

28. W. F. Stinespring, Positive Functions on C*-algebras, Proc. Am. Math. Soc. 6, 211 (1955). 

29. A. Jamiolkowski, Linear transformations which preserve trace and positive semidefiniteness 
of operators. Rep. Math. Phys. 3, 275 (1972). 



Many-body physics from a quantum information perspective 



47 



30. M. Horodecki, P. Horodecki, and R. Horodecki, Separability ofn-particle mixed states: nec- 
essary and sufficient conditions in terms of linear maps, Phys. Lett. A 283, 1 (2001). 

31. S. L. Woronowicz, Positive maps of low dimensional matrix algebras. Rep. Math. Phys. 10, 
165 (1976). 

32. K. Tanahashi and J. Tomiyama, Indecomposable positive maps in matrix algebras. Can. 
Math. Bull 31. 308 (1988). 

33. M. Horodecki and P. Horodecki, Reduction criterion of separability and limits for a class of 
distillation protocols, Phys. Rev. A 59, 4206 (2000). 

34. N. J. Cerf, C. Adami, and R. M. Gingrich, Reduction criterion for separability, Phys. Rev. A 
60, 898 (1999). 

35. H.-P. Breuer, Optimal Entanglement Criterion for Mixed States, Phys. Rev. Lett. 97, 080501 
(2006). 

36. W. HaR, A new criterion for indecomposability of positive maps, J. Phys. A 39, 14119 (2006). 

37. S. Woronowicz, Nonextendible positive maps. Comm. Math. Phys. 51, 243 (1976). 

38. M. Lewenstein, B. Kraus, J. 1. Cirac, and P. Horodecki, Optimization of entanglement wit- 
nesses, Phys. Rev. A 62, 052310 (2000). 

39. B. M. Terhal, A family of indecomposable positive linear maps based on entangled quantum 
states, Lin. Alg. Appl. 323, 61 (2001). 

40. C. H. Bennett, D. P. DiVincenzo, J. A. Smolin, and W. K. Wootters, Mixed-state entanglement 
and quantum error correction, Phys. Rev. A 54, 3824 (1996). 

41. V. Vedral, M. B. Plenio, M. A. Rippin, and P. L. Knight, Quantifying Entanglement, Phys. 
Rev. Lett 78, 2275 (1997). 

42. D. P. DiVincenzo, C. A. Fuchs, H. Mabuchi, J. A. Smolin, A. Thapliyal, and A. Uhhnann, in 
Proceedings of the first NASA International Conference on Quantum Computing and Quan- 
tum Communication (Springer- Verlag, 1998). 

43. T. Laustsen, F. Verstraete, and S. J. van Enk, Local vs. joint measurements for the entangle- 
ment of assistance. Quantum Inf. Comput. 3, 64 (2003). 

44. M. A. Nielsen, Conditions for a Class of Entanglement Transformations, Phys. Rev. Lett. 83, 
436(1999). 

45. G. Vidal, Entanglement monxttones, 1. Mod. Opt. 47, 355 (2000). 

46. D. Jonathan and M. B. Plenio, Minimal Conditions for Local Pure-State Entanglement Ma- 
nipulation, Phys. Rev. Lett. 83, 1455 (1999). 

47. M. Horodecki, A. Sen(De), and U. Sen, Dual entanglement measures based on no local 
cloning and no local deleting, Phys. Rev. A 70, 052326 (2004). 

48. M. Horodecki, Distillation and bound entanglement. Quantum Inf. Comput. 1, 3 (2001). 

49. M. B. Plenio and S. Virmani, An introduction to entanglement measure. Quant. Inf. Comp. 
7, 1 (2007). 

50. C. H. Bennett, H. J. Bernstein, S. Popescu, and B. Schumacher, Concentrating partial entan- 
glement by local operations, Phys. Rev. A 53, 2046 (1996). 

51. S. Hill and W. K. Wootters, Entanglement of a Pair of Quantum Bits, Phys. Rev. Lett. 78, 
5022 (1997). 

52. W. K. Wootters, Entanglement of Formation of an Arbitrary State of Two Qubits, Phys. Rev. 
Lett. 80, 2245 (1998). 

53. B. M. Terhal and K. G. H. VoUbrecht, Entanglement of Formation for Isotropic States, Phys. 
Rev. Lett. 85, 2625 (2000). 

54. K. G. H. Vollbrecht and R. F. Werner, Entanglement measures under symmetry, Phys. Rev. 
A 64, 062307 (2001). 

55. P. Rungta, V. Buzek, C. M. Caves, M. Hillery, and G. J. Milburn, Universal state inversion 
and concurrence in arbitrary dimensions, Phys. Rev. A 64, 042315 (2001). 

56. P. Rungta and C. M. Caves, Concurrence-based entanglement measures for isotropic states, 
Phys. Rev. A 67, 012307 (2003). 

57. L. Aolita and F. Mintert, Measuring Multipartite Concurrence with a Single Factorizable 
Observable, Phys. Rev. Lett. 97, 050501 (2006). 

58. S. P. Walbom, P. H. S. Ribero, L. Davidovich, F. Mintert, and A. Buchleitner, Experimental 
determination of entanglement with a single measurement. Nature 440, 1022 (2006). 



48 



R. Augusiak, F. M. Cucchietti, and M. Lewenstein 



59. G. Vidal and R. F. Werner, Computable measure of entanglement, Phys. Rev. A 65, 032314 
(2002). 

60. M. B. Plenio, Logarithmic Negativity: A Full Entanglement Monotone That is not Convex, 
Phys. Rev. Lett. 95, 090503 (2005). 

61. L. Bombelli, R. K. Koul, J. Lee, and R. D. Sorkin, Quantum source of entropy for black 
holes, Phys. Rev D 34, 373 (1986). 

62. M. Srednicki, Entropy and area, Phys. Rev. Lett. 71, 666 (1993). 

63. J. D. Bekenstein, Black Holes and Entropy, Phys. Rev. D 7, 2333 (1973). 

64. J. D. Bekenstein, Black holes and information theory, Contemp. Phys. 45, 31 (2004). 

65. S. W. Hawking, Black hole explosions?. Nature 248, 30 (1974). 

66. R. Bousso, The holographic principle. Rev. Mod. Phys. 74, 825 (2002). 

67. J. Eisert, M. Cramer, and M. B. Plenio, Area laws for the entanglement entropy - a review. 
Rev. Mod. Phys. 82, 277 (2010). 

68. P. Calabrese, J. Cardy, and B. Doyon, Special Issue: Entanglement entropy in extended quan- 
tum systems, J. Phys. A 42, 500301 (2009). 

69. E. Lubkin, Entropy of an n-system from its correlation with a k-reservoir, J. Math. Phys. 
19, 1028 (1978). 

70. S. Lloyd and H. Pagels, Complexity as thermodynamic depth, Ann. Phys. (N.Y.) 188, 186 
(1988). 

71. D. N. Page, Average Entropy of a Subsystem, Phys. Rev. Lett. 71, 1291 (1993). 

72. I. Bengtsson and K. Zyczkowski, Geometry of Quantum States (Cambridge University Press, 
2006). 

73. S. K. Foong and S. Kanno, Proof of a Page's Conjecture on the Average Entropy of a Sub- 
system, Phys. Rev Lett. 72, 1148 (1994). 

74. S. Sen, Average Entropy of a Quantum Subsystem, Phys. Rev. Lett. 77, 1 (1996). 

75. J. Sanchez-Ruiz, Simple proof of Pages conjecture on the average entropy of a subsystem, 
Phys. Rev E 52, 5653 (1995). 

76. M. B. Hastings, An area law for one-dimensional quantum system, J. Stat. Mech. Theory and 
Exp. p. P08024 (2007). 

77. E. H. Eisert and D. W. Robinson, The finite group velocity of quantum spin systems, Conmi. 
Math. Phys. 28, 251 (1972). 

78. L. Masanes, Area law for the entropy of low-energy states, Phys. Rev. A 80, 052104 (2009). 

79. W. Diir, L. Hartmaim, M. Hein, M. Lewenstein, and H.-J. Briegel, Entanglement in Spin 
Chains and Lattices with Long-Range Ising-Type Interactions, Phys. Rev. Lett. 94, 097203 

(2005) . 

80. J. Eisert and T. Osborne, General Entanglement Scaling Laws from Time Evolution, Phys. 
Rev Lett. 97, 150404 (2006). 

81. J. 1. Latorre and A. Riera, A short review on entanglement in quantum spin systems, J. Phys. 
A 42, 504002 (2009). 

82. G. Vidal, J. I. Latorre, E. Rico, and A. Kitaev, Entanglement in Quantum Critical Phenom- 
ena, Phys. Rev Lett. 90, 227902 (2003). 

83. B.-Q. Jin and V. E. Korepin, Quantum Spin Chain, Toeplitz Deteminants and the Ficher- 
Hartwig Conjecture, J. Stat. Phys. 116, 79 (2004). 

84. A. R. Its, B.-Q. Jin, and V. E. Korepin, Entanglement in the XY spin chain, J. Phys. A: Math. 
Gen. 38, 2975 (2005). 

85. J. P. Keating and F. Mezzadri, Entanglement in Quantum Spin Chains, Symmetry Classes of 
Random Matrices, and Conformal Field Theory, Phys. Rev. Lett. 94, 050501 (2005). 

86. J. Eisert and M. Cramer, Single-copy entanglement in critical quantum spin chains, Phys. 
Rev A 72, 042112 (2005). 

87. P. Calabrese and J. Cardy, Entanglement entropy and conformal field theory, J. Phys. A 42, 
504005 (2009). 

88. M. M. Wolf, Violation of the Entropic Area Law for Fermions, Phys. Rev. Lett. 96, 010404 

(2006) . 

89. D. Gioev and I. Klich, Entanglement Entropy of Fermions in Any Dimension and the Widom 
Conjecture, Phys. Rev. Lett. 96, 100503 (2006). 



Many-body physics from a quantum information perspective 



49 



90. S. Farkas and Z. Zimboras, The von Neumann entropy asymptotics in multidimensional 
fermionic systems, J. Math. Phys. 48, 102110 (2007). 

91. M. B. Hastings, Locality in Quantum and Markov Dynamics on Lattices and Networks, Phys. 
Rev. Lett. 93, 140402 (2004). 

92. S. Boyd and L. Vanderberghe, Convex Optimization (Cambridge University Press, 2004). 

93. B. Groisman, S. Popescu, and A. Winter, Quantum, classical, and total amount of correla- 
tions in a quantum state, Phys. Rev. A 72, 032317 (2005). 

94. M. M. Wolf, F. Verstraete, M. B. Hastings, and J. I. Cirac, Area Laws in Quantum Systems: 
Mutual Information and Correlations, Phys. Rev. Lett. 100, 070502 (2008). 

95. F. Verstraete, M. Popp, and J. L Cirac, Entanglement versus Correlations in Spin Systems, 
Phys. Rev. Lett. 92, 027901 (2004). 

96. G. Vidal, Efficient Classical Simulation of Slightly Entangled Quantum Computations, Phys. 
Rev. Lett. 91, 147902 (2003). 

97. D. Perez-Garci'a, F. Verstraete, M. M. Wolf, and J. I. Ciarac, Matrix product state represen- 
tation. Quantum Inf. Comput. 7, 401 (2007). 

98. F. Verstraete and J. 1. Cirac, Renormalization algorithms for Quantum-Many Body Systems 
in two and higher dimensions, cond-mat/0407066 (2004). 

99. N. Schuch, M. M. Wolf, F. Verstraete, and J. I. Cirac, Computational Complexity of Projected 
Entangled Pair States, Phys. Rev. Lett. 98, 140506 (2007). 

100. I. Affleck, T. Kennedy, E. H. Lieb, and H. Tasaki, Rigorous results on valence-bond ground 
states in antiferromagnets, Phys. Rev. Lett. 59, 799 (1987). 

101. 1. Affleck, T. Kennedy, E. H. Lieb, and H. Tasaki, Valence bond ground states in isotropic 
quantum antiferromagnets, Commun. Math. Phys. 115, 477 (1988). 

102. C. K. Majumdar and D. K. Ghosh, On Next-Nearest-Neighbor Interaction in Linear Chain. 
I, J. Math. Phys. 10, 1388 (1969). 

103. R. J. Glauber, Time-Dependent Statistics of the Ising Model, J. Math. Phys 4, 294 (1963). 

104. U. Deker and F. Haake, Renormalization Group Transformation for the Master Equation of 
a Kinetic Ising Chain, Z. Phys. B 35, 281 (1979). 

105. J. C. Kimball, The kinetic Ising model: Exact susceptibilities of two simple examples, J. Stat. 
Phys. 21, 289 (1979). 

106. F. Haake and K. Thol, Universality Classes for One Dimensional Kinetic Ising Models, Z. 
Phys. 6 40,219(1980). 

107. B. U. Felderhof, Spin relaxation of the Ising chain. Rep. Math. Phys. 1, 215 (1971). 

108. E. D. Siggia, Pseudospin formulation of kinetic Ising models, Phys. Rev. B 16, 2319 (1977). 

109. S. R Heims, Master Equation for Ising Model, Phys. Rev. 138, A587 (1965). 

1 10. K. Kawasaki, in Phase Transition and Critical Phenomena, vol. 2, ed. by C. Domb and M. S. 
Green (Academic Press, 1972)., chap. 11, pp. 443-501 

111. R. Augusiak, F. M. Cucchietti, F. Haake, and M. Lewenstein, Quantum kinetic Ising models. 
New J. Phys. 12, 025021 (2010). 

112. H. J. Hilhorst. M. Suzuki, and B. U. Felderhof, Kinetics of the stochastic Ising chain in a 
two-flip model, Physica 60, 199 (1972). 

113. P. Jordan and E. Wigner, fiber das Paulische Aequivalenzverbot, Z. Phys. 47, 631 (1928). 

114. N. N. Bogoliubov, On a new method in the theory of superconductivity, Nuovo Cimento 7, 
794(1958). 

115. J. G. Valatin, Comments on the theory of superconductivity, Nuovo Cimento 7, 843 (1958). 

116. L. Amico, R. Fazio, A. Osterloh, and V. Vedral, Entanglement in many-body systems. Rev. 
Mod. Phys. 80, 517 (2008). 

117. P. Zanardi and N. Paunkovic, Ground state overlap and quantum phase transitions, Phys. 
Rev. E 74, 031123 (2006). 

118. H. T. Quan, Z. Song, X. F. Liu, P. Zanardi, and C. P. Sun, Decay of Loschmidt Echo Enhanced 
by Quantum Criticality, Phys. Rev. Lett. 96, 140604 (2006). 

119. C. Zhang, S. Tewari, R. Lutchyn, and S. D. Sarma, px+ipy Superfluid from s-Wave Interac- 
tions of Fermionic Cold Atoms, Phys. Rev. Lett. 101, 160401 (2008). 



50 



R. Augusiak, F. M. Cucchietti, and M. Lewenstein 



120. J. Zhang, F. M. Cucchietti, C. M. Chandrashekar, M. Laforest, C. A. Ryan, M. Ditty, A. Hub- 
bard, J. K. Gamble, and R. Laflamme, Direct observation of quantum criticality in Ising spin 
chains, Phys. Rev. A 79, 012305 (2009). 

121. H. Li and F. D. M. Haldane, Entanglement Spectrum as a Generalization of Entanglement 
Entropy: Identification of Topological Order in Non-Abelian Fractional Quantum Hall Effect 
States, Phys. Rev. Lett. 101, 010504 (2008). 

122. P. Calabrese and A. Lefevre, Entanglement spectrum in one-dimensional systems, Phys. Rev. 
A 78, 032329 (2008). 

123. G. Vidal, Entanglement Renormalization, Phys. Rev. Lett. 99, 220405 (2007). 

124. S. R. Clark and D. Jaksch, Dynamics of the superfluid to Mott-insulator transition in one 
dimension, Phys. Rev. A 70, 043612 (2004). 

125. C. V. Kraus, N. Schuch, F. Verstraete, and J. I. Cirac, Fermionic projected entangled pair 
states, Phys. Rev. A 81, 052338 (2010). 

126. P. Corboz and G. Vidal, Fermionic multiscale entanglement renormalization ansatz, Phys. 
Rev. B 80, 165129 (2009). 

127. P. Corboz, G. Evenbly, F. Verstraete, and G. Vidal, Simulation of interacting fermions with 
entanglement renormalization, Phys. Rev. A 81, 010303 (2010). 

128. T. Barthel, C. Pineda, and J. Eisert, Contraction of fermionic operator circuits and the simu- 
lation of strongly correlated fermions, Phys. Rev. A 80, 042333 (2009). 

129. P. Corboz, R. Orus, B. Bauer, and G. Vidal, Simulation of strongly correlated fermions in two 
spatial dimensions with fermionic projected entangled-pair states, Phys. Rev. B 81, 165104 
(2010). 

130. C. Pineda, T. Barthel, and J. Eisert, Unitary circuits for strongly correlated fermions, Phys. 
Rev. A 81, 050303 (2010). 

131. M. Troyer and U.-J. Wiese, Computational Complexity and Fundamental Limitations to 
Fermionic Quantum Monte Carlo Simulations, Phys. Rev. Lett. 94, 170201 (2005). 

132. K. Kim, M.-S. Chang, S. KorenbUt, R. Islam, E. E. Edwards, J. K. Freericks, G.-D. Lin, L.- 
M. Duan, and C. Monroe, Quantum simulation of frustrated Ising spins with trapped ions. 
Nature 465, 590 (2010). 

133. R. Jordens, L. Tarruell, D. Greif, T. Uehlinger, N. Strohmaier, H. Moritz, T. Esslinger, 

L. De Leo, C. Kollath, A. Georges, V. Scarola, L. PoUet, E. Burovski, E. Kozik, and 
M. Troyer, Quantitative Determination of Temperature in the Approach to Magnetic Order 
of Ultracold Fermions in an Optical Lattice, Phys. Rev. Lett. 104, 180401 (2010). 

134. K. Temme, M. M. Wolf, and F. Verstraete, Stochastic exclusion processes versus coherent 
transport, e-print arXiv:0912.0858 (2009). 

135. F. Verstraete, M. M. Wolf, and J. I. Cirac, Quantum computation and quantum-state engi- 
neering driven by dissipation, Nat. Phys. 5, 633 (2009). 

136. B. Kraus, H. P. Biichler, S. Diehl, A. Kantian, A. Micheli, and P. ZoUer, Preparation of 
entangled states by quantum Markov processes, Phys. Rev. A 78, 042307 (2008). 

137. S. Diehl, A. Micheli, A. Kantian, B. Kraus, H. P. Biichler, and P. Zoller, Quantum states and 
phases in driven open quantum systems with cold atoms, Nat. Phys. 4, 878 (2008). 



